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1*.  AMTRACT  (Gamma  on  rouono  4  lawny  JAW  Moootf  Ry  RRRCR  nomOor) 

>This  report  describes  theoretical  models  of  laser/plasma  ex perJmentSy currently  under  way 
at  the  Naval  Research  Laboratory  (NRL)  <*in  which  hot  laser- irradiated  target  debris 
expands  into  a  few  Torr  background  gas.  The  most  interesting  characteristics  of  these 
experiments  are  the  formation  and  evolution  of  a  spherical  blast  wave,  and  in  some  cases, 
the  appearance  of  one  or  more  distinctive  protuberances  or  ^aneurysms^  on  the  expanding 
wave.  Many  possible  mechanisms  for  the  formation  of  these  structures  have  been  sug¬ 
gested,  ranging  from  instabilities  in  the  blast  wave  itself  to  the  effects  of  jets  of 
material  ejected  from  the  target.  Many  of  these  mechanisms  are  of  interest  to  the 
problem^  of  early-time  structuring  in  high  altitude  nuclear  bursts.  ■Tnr"th1s  report  we 
discuss^a  series  of  magnetohydrodynamic  models  that  accurately  describe  this  phenomenon, 
and  identify  the  processes  that  mold  the  blast  wave  and  the  aneurysm.  We  sfanw  t hatpin 
most  cases  the  aneurysm  can  be  explained  by  a  channel  heating  mechanism  due  to  inverse 
bremsstrahlung  absorption  in  the  laser  channel.  A  simple  model  is  derived  that  clearly 
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shows  this  effect  and  matches  the  detailed  two-dimensional  models  as  well  as  the  experi¬ 
mental  results.  The  effects  of  laser-target  interaction,  thermal  conduction,  radiation 
cooling,  multiple  materials,  self-generated  magnetic  fields,  and  ambient  magnetic  fields 
are  included  and  discussed. (  ! A  i 
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We  have  benefited  from  many  fruitful  discussions  with  H.  Frese 
and  J.  Ruff  of  MRC,  and  with  the  Naval  Research  Laboratory  (NRL)  experi¬ 
mentalists,  especially  R.  H.  Ripin  and  J.  A.  Stamper  (who  kindly  provided 
us  with  the  shadowgraphs). 
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SECTION  1 


INTRODUCTION 

From  1983  through  the  present,  the  NRL  laser  group  has  performed 
a  large  number  of  experiments  that  are  of  interest  to  the  HANE  community 
with  the  Pharos  II  and  Pharos  III  Nd-laser  and  target  facilities.  The 
maximum  laser  energy  for  the  earlier  Pharos  II  experiments  was  on  the  order 
of  100  J,  whereas  the  energy  in  experiments  of  interest  to  us  with  the  new 
Pharos  III  facility  ranged  from  5  to  300  0  in  each  of  up  to  three  beams. 

In  most  of  these  experiments,  a  single  beam  from  the  Nd-glass  laser  A  = 
1.054  um,  tr  ,  -4  ns)  is  focused  onto  a  solid  planar  aluminum  or 

polystyrene  (CH)  target  that  is  immersed  in  an  ambient  medium  (such  as 
nitrogen)  at  a  pressure  of  between  a  few  mTorr  and  10  Torr.  The  laser 
focal  spot  diameter  is  typically  between  0.2  and  1.0  mm,  and  the  targets 
are  thin,  usually  between  4  and  25  urn  thick.  In  the  subset  of  those 
experiments  designated  as  being  in  the  "collisional"  or  "high-pressure" 
regime,  the  ambient  pressure  is  above  1  Torr.  The  laser,  with  irradiance 
between  1012  and  1014  Watt/cm2,  ablates  target  material  and  imparts  to  the 
ions  a  controllable  velocity  distribution  (0.1  <  au/u  <  1.0)  around  a  peak 
of  a  few  107  cm/sec  (Grun,  Stel 1 ingwerf ,  and  Ripin,  1986,  Ref.  1).  In  the 
"high-pressure"  regime,  expanding  target  plasma  couples  to  the  ambient 
material  so  that,  on  time  scales  long  compared  to  the  laser  pulse  length,  a 
thin,  almost  spherical  blast  front  is  observed  (Ripin,  et  al.,  1986,  Ref. 
?).  Because  the  debris  expands  in  the  direction  of  the  laser  beam,  the 
blast  front  is  not  centered  about  the  target  but  rather  exDands  prefer¬ 
entially  toward  the  source  of  the  laser  light. 

A  complete  tabulation  of  the  NRL  laser  shots  through  early  1985 
is  available  in  Reference  3  (Stel 1 ingwerf ,  1987),  including  all  shot  param¬ 
eters  and  blast  wave  configurations. 

A  variety  of  diagnostic  tools  were  used  to  investigate  the  laser- 
target-ambient  interactions  (Stamper,  et  al.,  1986,  Ref.  4).  Although  the 
’"adius  of  the  expanding  shell  scales  as  a  classical  blast  wave  (Ripin, 
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Stamper,  McLean,  1984,  Ref.  5),  a  number  of  non-blast-wave  features  are 
evident  from  the  diagnostics.  We  have  modeled  the  laser-target-ambient 
interactions  with  the  two-dimensional  MHD  computer  code,  MACH2  (Frese,  1987 
and  Peterkin,  et  al.,  1986,  Refs.  6  and  7),  and  are  able  to  simulate  the 
formation  and  time-development  of  the  nearly-spherical  blast  wave 
(Stel 1 ingwerf ,  Peterkin,  Longmire,  1986,  Ref.  8). 

The  most  prominent  non-blast-wave  feature  of  the  expanding  bubble 
seen  in  diagnostics  of  the  "high  pressure"  experiments  is  the  aneurysm.  It 
appears  as  an  outward-pointing  protrusion  attached  to  the  otherwise  spheri¬ 
cal  blast  wave.  In  shadowgraphs  from  experiments  done  with  the  Pharos  II 
facility,  the  aneurysm  appeared  at  various  locations  within  an  imaginary 
cone  of  about  30  degree  half  angle  that  was  centered  about  the  target 
normal.  The  laser  beam  was  always  at  normal  incidence  to  the  target. 
Because  the  aneurysm  appeared  to  be  aligned  with  neither  the  the  laser  beam 
nor  the  target  normal,  the  mechanism  responsible  for  its  formation  was  not 
obvious. 


To  explain  the  aneurysm,  a  number  of  interesting  ideas  were  put 
forth.  These  include: 

1.  non-uniform  laser  intensity  across  the  laser  spot, 

2.  a  jet  of  material  from  the  target  that  catches  up  with  the 
leading  blast  front  (as  suggested  by  C.  Longmire), 

3.  instability  of  blast  wave  induced  by  thermally  generated 
magnetic  fields 

4.  laser  heating  of  the  plasma  within  the  laser  path. 

We  modeled  the  first  of  these  processes  with  MACH2,  and  concluded 
that  it  was  not  the  cause  of  the  aneurysm  (Stel 1 ingwerf ,  Peterkin, 

Longmire,  1986,  Ref.  8).  Subsequently,  experiments  were  done  at  NRL  with 
the  Pharos  III  facility  that  employed  an  ISI-smoothed  laser  pulse 
(Lehmberg,  Obenschain,  1983,  Ref.  9).  The  aneurysm  persisted,  thus  ruling 
out  a  spatially  non-uniform  pulse  as  the  driving  mechanism.  As  for  the 
second  idea,  the  material  that  is  ablated  from  the  target  at  late  times 
travels  through  a  region  from  which  ambient  gas  has  been  swept  away  by  the 
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leading  blast  front.  The  leading  front  may  be  overtaken  by  later  debris. 
This  may  explain  aneurysms  that  appear  normal  to  the  target  when  it  is 
illuminated  by  a  low  energy  laser  pulse  with  non-normal  incidence.  Such  is 
seen  in  one  of  the  three  NRL  experiments  done  at  very  low  energy  with  the 
Pharos  III  facility.  (Experiments  86-582,  -583,  and  -586  were  low  energy 
shots;  86-583  shows  the  feature.)  We  have  modeled  thermally  generated 
magnetic  fields,  the  third  from  this  list,  with  MACH2  (Peterkin, 

Stel 1 ingwerf,  1987,  Ref.  10).  In  a  model  with  parameters  from  NRL  shot 
85-544  (109  J  with  5  Torr  ambient  nitrogen),  we  calculate  that  less  than 
one  percent  of  the  laser  energy  is  converted  to  magnetic  field  energy  near 
the  target,  and  this  field  is  convected  with  the  expanding  bubble.  Our 
simulations  show  no  disruption  of  the  shell  from  this  magnetic  pertur¬ 
bation.  Experimental  work  is  under  way  to  investigate  the  magnitude  and 
location  of  thermally  generated  magnetic  fields.  It  has  been  suggested 
(Suillani,  1985,  Ref.  11)  that  such  fields  may  explain  the  turbulent- 
looking  structure  that  seems  to  exist  in  front  of  the  blast  wave  as  seen  in 
shadowgraphs  of  some  experiments  with  low  energy  laser  pulses. 

Clearly,  the  second  and  third  ideas  from  the  above  list  deserve 
further  consideration  in  attempts  to  explain  deviations  from  classical 
blast  wave  behavior  seen  in  laser-target-ambient  experiments.  In  this 
paper,  however,  shall  limit  our  discussions  to  the  laser  heated  channel. 

We  illustrate  that  if  the  laser  energy  is  sufficient,  an  aneurysm  will  form 
centered  about  the  laser  axis.  This  idea  has  been  investigated  previously 
(Lyon,  1985,  Ref.  12)  with  a  null  result.  Our  two-dimensional  modeling 
shows,  however,  that  laser  heating  of  a  narrow  channel  does  cause  an 
aneurysm  when  the  laser  energy  is  at  least  a  few  tens  of  Joules  and  when 
the  ambient  pressure  is  above  a  Torr  ( Stel 1 ingwerf ,  Peterkin,  Longmire, 
1986,  Ref.  8). 

Although  an  aneurysm  is  always  present  centered  about  the  laser 
axis  when  the  laser  energy  is  above  a  threshold  on  the  order  of  10  J  in 
experiments  with  the  new  Pharos  III  facility,  small  aneurysm-like  Dro- 
trusions  are  often  visible  off-axis  as  well.  It  seems  likely  that  there  is 
nore  than  a  single  mechanism  responsible  for  the  production  of  such 
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features  on  the  shell.  However,  the  experimental  reproducibi 1 ity  of  the 
on-axis  aneurysm,  and  our  ability  to  simulate  it  by  modeling  a  laser-heated 
channel  in  the  real  experimental  geometry,  suggest  that  the  laser  heating 
of  the  photoionized  plasma  in  the  channel  is  the  driving  mechanism  for  the 
on- ax is  aneurysm. 

Structure  similar  to  the  aneurysm  is  seen  in  photographs  of  the 
blast  from  nuclear  devices  that  are  suspended  from  tethered  balloons  or 
support  towers  and  stabilized  by  guy-wires  (Jones,  Childers,  Kodis,  1987, 
Ref.  13).  Called  "cable  spikes",  these  conical-shaped  protrusions  emanate 
from  the  expanding  blast  front  along  the  guy-wires.  Cable  spikes  are 
thought  to  be  caused  by  absorption  of  radiation  by  the  wires  and  the 
development  of  a  heated,  low  density  channel  (similar  to  the  present  case). 

In  this  report  the  aneurysm  is  discussed  primarily  from  the  view¬ 
point  of  the  channel  heating  mechanism,  but  other  physical  processes  are 
addressed  as  well.  Section  2  presents  motivation,  background,  and  results 
of  several  series  of  two-dimensional  MHD  simulations  that  demonstrate 
aneurysm  formation  and  its  dependence  on  various  parameters.  In  Section  3 
a  simple  channel  heating  model  is  developed  that  explains  the  detailed  MHD 
results.  Section  4  discusses  the  method  of  computing  sel f-oenerated 
magnetic  fields  used  in  the  code,  while  Section  5  shows  a  their  effects  in 
sample  computations.  Finally,  Section  6  shows  a  preliminary  computation  of 
a  case  with  both  self-generated  fields  and  an  ambient  magnetic  field 
parallel  to  the  laser  path. 

A  summary  of  the  experimental  and  theoretical  results  of  this 
investigation  have  been  submitted  to  The  Physics  of  Fluids  and  should 
appear  in  late  1988  or  early  1989. 
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SECTION  2 


THE  ANEURYSM 

2.1  THE  CHANNEL-HEATING  MECHANISM. 

A  simple  calculation  illustrates  how  the  temperature  of  the  gas 
within  the  path  of  a  high-irradiance  laser  can  rise  from  room  temperature 
to  many  electron  Volts.  Let  us  assume  that  X-rays  emanating  from  the 
1 aser-irradiated  target  singly  ionize  a  fraction  of  the  constituent  atoms 
of  the  gas  nearby  during  the  early  stages  of  the  laser  pulse.  Subsequent 
laser  light  is  absorbed  by  the  ionized  atoms  near  the  target,  raising  the 
local  temperature  substantially.  This  hot  plasma  close  to  the  target  then 
itself  radiates,  ionizing  ambient  gas  nearby  in  the  process.  This  process 
cascades  outward  from  the  focal  spot  of  the  laser.  We  approximate  this 
complicated  radiation  transport  process  by  simply  assuming  the  ambient  gas 
is  preionized  in  a  manner  that  falls  as  the  square  of  the  distance  from  the 
focal  spot.  Such  1/r2  behavior  follows  from  the  SDherical  divergence  of 
the  laser  beam.  To  be  orecise,  our  model  for  oreioni zation  ;s  used  as  a 
’ower  bound  on  the  ionization  fraction  of  the  ambient  medium,  ^he  'oni z- 
ation  fraction,  c,  decreases  with  the  square  of  the  distance  from  a  ooint 
at  the  focal  soot  of  the  laser  and  is  normalized  to  have  the  value,  ;0,  at 
a  distance  one  cm  from  that  point.  We  tried  various  values  of  Co  and 
conclude  that  a  value  of  five  percent  gives  results  most  consistent  with 
experiment.  This  value,  however,  is  thought  to  be  higher  than  that  of 
•"ea!  ity. 


In  a  distance,  dr,  some  fraction  of  the  laser  light  is  absorbed 
by  the  free  electrons  via  inverse  bremsstrahlung.  The  fractional  power 
absorption  (Johnston  and  Dawson,  1973,  Ref.  14),  f,  in  mks  units  with 
temoerature  in  eV,  is: 
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where  a  is  the  constant  1.0779  10“5,  b  is  taken  to  be  0.5,  c*  is  the  laraer 
of  the  fractional  ionization,  c,  and  unity.  In  A  is  a  Coulomb  loqarithm  (of 
order  unity),  X  is  the  laser  wavelength  in  meters,  T  is  the  temoerature  in 
eV,  x  is  the  dimensionless  ratio  of  the  electron  number  density  to  the 
critical  number  density,  and  dr  is  an  increment  in  the  direction  of  laser 
propagation.  The  laser  light  that  arrives  at  the  critical  surface  is 
absorbed  there.  If  the  integrand  in  Equation  (1)  is  small  and  x  €  1  (that 
is  the  case  in  low  density  regions),  the  absorption  coefficient  is: 


be'  In  A  x' 

a - T7T72 — 


X“  T 

which  can  be  written  as 


«  u  >2  tr  p-  In  A 

sox  - 1  "m — 

Ac  ji/i 


(2) 


(3) 


where  S  is  3.133  1018,  p  is  the  mass  density,  and  A  is  the  atomic  weiqht  of 

the  plasma.  If  a  fraction,  <  ,  of  the  incident  laser  power,  P  ,  that 

v  1 

passes  through  a  disk  of  radius  r  and  density  p  is  absorbed  per  unit  length 
by  that  disk,  the  temoerature  change  is: 


dT 

HT  = 


K  ?_ 
V  0 


p  it  r  c, 


(4) 


where  c^  is  the  specific  heat  of  the  material.  If  we  assume  an  ideal  gas 
eauation  of  state,  the  specific  heat  is 


cy  «  1.44  10 


3  (i+c) 


Joule 

eV-kg 


(5) 
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independent  of  the  temperature,  and  we  can  estimate  the  size  of  the  temper¬ 
ature  rise  for  some  particular  values  of  the  parameters.  Equation  5  does 
not  take  into  consideration  the  energy  that  ionizes  rather  than  heats  the 
plasma.  The  ionization  energy  is  in  general  not  negligible  (see  Section  3, 
Figure  20).  Combining  Equations  (3),  (4),  and  (5),  we  have 


dT 

■3T 


y  b  x  c*  c  p  In  a  Pj 
r2  A  T3/2  (1  +  0  " 


(6) 


where  y  =  6.925  109.  If  c',  c,  p,  In  A,  and  r  are  approximated  as  time 
independent  quantities,  we  can  integrate  Equation  (6).  The  time  integral 
of  the  incident  laser  power  is  the  incident  laser  energy,  E0.  We  then  have 
for  the  final  electron  temperature,  T^,  after  a  time  interval,  t 


rT  5/2  A  5  f  Y  ^  C'  ^  p  ln  A  b  E0  n2/5 

L 1  o  7  ^ - ? - ; - JJ 

0  2  r2  A  (1  +  c) 


(7) 


where  T0  is  the  initial  temperature. 

Let  us  try  to  estimate  the  temperature  rise  in  the  laser  channel 
for  NRL  experiment  85-544.  Oark  field  shadowgraphs  from  this  experiment 
show  a  well-formed  aneurysm  centered  along  the  laser  axis  as  illustrated  in 
Figure  1.  A  109  J  laser  beam  propagates  through  room  temperature  Nitrogen 
gas  at  5  Torr  pressure.  This  corresponds  to  a  mass  density  of  8.2 
10-6  gm/cc,  which  if  five  percent  ionized  (five  percent  of  the  ions  are 
singly  ionized)  at  a  distance  of  1  cm  above  that  focal  spot  of  the  laser 
produces  an  electron  number  density  of  8.77  1015/cc.  The  laser  is  focused 
to  a  circular  spot  approximately  100  ym  in  radius,  and  the  beam  converges 
to  this  spot  with  an  angle  of  2.86  degrees  from  the  target  normal  (as 
illustrated  in  Figure  2).  At  a  distance  one  cm  from  the  focal  spot,  the 
radius  of  the  laser  disk  is  almost  600  um  and  the  cross  sectional  area  of 
the  ambient  gas  that  is  pierced  by  the  laser  beam  there  is  a  larger  than 
the  laser  focal  spot  by  a  factor  of  36.  During  the  first  half  of  the 
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Figure  1.  Dual  time  shadowgram  (71  ns  and  140  ns)  of  WL  85-544. 
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Figure  2a.  Outline  of  the  calculation  mesh  for  MACH2  simulation,  showing 
dimensions.  The  target  region  is  enlarged  in  Figure  2b. 


Figure  2b.  Outline  of  the  target  region  for  MACH2  simulation,  showing 
dimensions. 
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pulse,  54.5  J  of  energy  is  incident  on  any  slice  of  the  ambient  gas  if 
f  «  1.  From  Equation  (7),  we  can  estimate  as  a  lower  limit  that  the 
temperature  of  a  plasma  slice  600  ym  in  radius  will  rise  from  room  tempera¬ 
ture  to  be  on  the  order  of  1.33  eV  if  an  ideal  gas  specific  heat  is  assumed 
and  if  the  ionization  and  mass  density  are  assumed  to  be  almost  constant. 

At  this  density,  nitrogen  is  fully  ionized  at  100  eV.  If  we  assume  full 
ionization  and  calculate  In  A  to  be  3.3,  we  find  an  upper  limit  of  108  eV 
on  the  temperature  1  cm  from  the  target.  Clearly,  a  self-consistant  model 
that  takes  the  ionization  into  account  is  required  to  accurately  estimate 
the  temperature  rise  in  the  channel.  Because  of  the  strong  focusing  of  the 
laser  beam,  laser  energy  density,  and  hence  the  change  in  temperature  of 
the  gas,  are  expected  to  vary  essentially  with  the  inverse  of  the  cross- 
sectional  area  of  the  beam.  If  we  use  our  simple  model  to  estimate  the 
rise  in  temperature  of  the  gas  at  the  focal  spot  of  the  laser--where  the 
intensity  is  36  times  as  large  as  it  is  1  cm  away- -we  obtain  roughly  3  KeV 
as  an  upper  limit  on  the  temperature  at  the  target. 

Our  estimate  for  the  fraction  of  the  incident  power  absorbed  1  cm 
from  the  laser  focal  spot  is  rough  because  as  the  temperature  rises  a  few 
tens  of  eV,  the  fractional  ionization  does  change  appreciably  (the  average 
nitrogen  atom  in  a  5  Torr,  5  eV  plasma  is  twice  ionized).  As  the  tempera¬ 
ture  rises  so  does  the  fractional  ionization,  but  as  the  heated  plasma 
expands  its  density  falls.  Because  the  fraction  of  the  laser  energy  that 
is  absorbed  by  the  medium  increases  as  the  fractional  ionization  rises,  but 
decreases  as  the  temperature  rises  and  the  density  falls,  a  linear  analysis 
is  not  sufficient  to  determine  the  exact  nature  of  the  laser  channel. 
Clearly,  though,  additional  heating  of  a  fully  ionized  material  will  serve 
to  decrease  the  fractional  absorption  as  the  temperature  increases. 

These  effects  are  included  in  a  more  comprehensive  model  in 
Section  3  of  this  report.  The  present  results  clearly  indicate  channel 
heating  will  be  an  important  process,  and  justify  detailed  numerical  models 
to  verify  this  conclusion. 
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2.2 


TWO-DIMENSIONAL  MODELS. 


We  have  performed  2  1/2  dimensional  computer  simulations  of  the 
NRL  laser-target-ambient  system  with  MACH2,  an  MHD  code  which  is  described 
in  detail  elsewhere  (Frese,  1987,  and  Peterkin,  et  al.,  1986,  Refs.  6  and 
7).  Briefly,  MACH2  is  a  sinqle  fluid,  single  temperature  MHD  code  that 
carries  all  three  spatial  components  of  vectors,  but  allows  no  quantity  to 
depend  on  the  azimuthal  coordinate  when  it  models  cyl indrically-symmetric 
problems.  Accurate  equations  of  state  that  assume  local  thermodynamic 
equilibrium  are  used  by  accessing  SESAME  tables  at  Los  Alamos  National 
Laboratory  (LANL)  during  each  computational  cycle.  The  ions  and  electrons 
are  assumed  to  be  in  equilibrium  and  hence  are  character! zed  by  a  single 
temperature.  The  net  charge  in  any  cell  is  zero.  The  code  solves  the 
following  equations  in  MKSA  units  in  a  fractional  time-split  manner: 

Magnetic  Induction 


if  ■ » « (» *  *) 


I  7  «  l(nT-5-)  eir  (3  * 


Fluid  Momentum 


3v 


p  7 r  =  '  p  v 


Specific  Internal  Energy 


|y=-pv*7e  +  (-PgJ1  +  aJ1)  v.  v.  +  miQ  J2 

"  eTn — i  'rT.T  ^  *  7  p  +  7  *  (|<7T)  “  3pxT4  +  f  ?1 
'  e  i ' 


(3) 


J  V1  +  Vj  |-Pg^  +  +  i-  (B^B 1  -  J  B2 


9) 


(10) 
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Mass  Continuity 


ff  =  -V  *  (pv)  (11) 

In  the  above,  B  is  the  magnetic  induction,  yD  the  permeability  of  free 
space,  v  the  fluid  velocity,  n  the  electric  diffusivity  (electric  resis¬ 
tivity  divided  by  u0),  ne  an(*  ni  the  electron  and  ion  number  densities,  P 
the  thermal  pressure,  p  the  mass  density,  aJ1  the  viscous  stresses,  e  the 

plasma  specific  internal  energy,  T  the  temperature,  J  the  current  density 
* 

((l/u0)  curl  8),  <  the  thermal  conductivity,  e  the  magnitude  of  the  charge 
of  an  electron,  gJ1  the  spatial  metric,  x  the  Planck  mean  opacity  (cross- 
section  per  unit  mass),  a  Stephan's  constant,  the  instantaneous  laser 
power  per  unit  volume  at  infinity,  and  f  the  fractional  absorption  given  by 
Equation  (1). 

The  code  is  axisymmetric  with  the  axis  of  symmetry  appearing  as  a 
vertical  line  to  the  left  in  computer-generated  pictures.  Rotation  about 
this  axis  is  assumed.  The  MACH2  simulations  discussed  herein  take  place  on 
numerically-generated  Eluerian  meshes.  The  difference  equations  are  of  the 
finite  volume  type  which  is  well  suited  to  non-rectangul ar  meshes  required 
for  problems  with  complex  shapes. 

Results  from  four  MACH2  calculations  of  a  focused  laser  pulse 
interacting  with  the  partially  ionized  gas  are  illustrated  in  Figure  3.  In 
this  figure  we  plot  the  temperature  on  a  log  scale  as  a  function  of 
distance  from  the  laser  focal  spot  through  a  slice  down  the  laser  channel. 
To  simplify  our  computations,  we  modeled  only  the  ambient  gas  (without  a 
solid  target)  in  each  of  the  simulations.  The  sole  difference  among  these 
simulations  is  the  degree  of  preionization.  Recall  that  our  model  for  pre¬ 
ionization  decreases  with  the  square  of  the  distance  from  that  focal  spot. 
Four  values  for  the  initial  fractional  ionization  at  1  cm  from  the  focal 
spot  were  chosen  (one  percent,  three  percent,  five  percent,  and  ten  per¬ 
cent).  The  temperature  rises  to  approximately  600  eV  at  the  laser  focal 
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Figure  3.  Temperature  at  4  ns  (time  of  peak  laser  power)  versus  distance  from 
focal  spot  for  a  5  Torr  preionized  plasma  (with  4  distinct  values) 
without  a  target. 


soot  (where  the  laser  intensity  is  6  1013  Watt/cm2)  in  each  of  these  simu¬ 
lations.  The  front  of  the  temperature  gradient  is  further  from  the  focal 
spot  as  the  amount  of  preionization  is  increased.  For  example,  54.5  J  from 
a  laser  raises  the  temperature  in  the  channel  1  cm  from  the  focal  spot  to 
almost  1  eV  when  the  initial  ionization  fraction  is  five  percent  there  (not 
inconsistent  with  the  previous  analysis),  but  to  only  0.25  eV  when  the 
initial  ionization  fraction  is  one  percent.  The  front  of  the  hot  channel 
extends  further  from  the  focal  spot  as  the  degree  of  preionization 
increases.  By  varying  the  preionization  in  simulations  with  a  tarqet  in 
place,  we  have  found  that  it  is  necessary  for  the  hot  channel  to  extend 
beyond  the  blast  wave  at  early  times  if  an  aneurysm  is  to  form. 

The  physics  is  further  complicated  when  a  dense  target  is  present 
to  absorb  the  1 arqe  fraction  of  the  laser  energy  that  is  not  absorbed  in 
the  ambient  gas.  Fiqures  4a-c  illustrate  a  simulation  of  MRL  shot  35-5*4, 
complete  with  a  tarqet  at  the  focal  SDot  of  the  laser,  where  we  have  ore- 
ionized  the  ambient  nitrogen  to  a  level  of  five  oercent  1  cm  from  the 
tarqet  surface.  As  seen  in  the  10  ns  picture,  Fiqure  *a,  the  hot  shell 
within  the  channel  (where  the  temperature  is  0.45  eV)  extends  1.53  cm  from 
the  target.  The  exoan ding  bubble  is  substantially  behind  at  this  time--a 
iistance  of  0.5  cm  from  the  taraet  surface.  The  tarqet  surface  is  not. 

(510  eV)  at  10  ns,  but  cools  to  60  eV  at  71  ns.  Figures  *b  (at  30  ns)  and 
4c  (71  ns)  illustrate  that  the  shell  gets  closer  to  the  temperature  front 
as  time  progresses.  However,  the  large  separation  at  early  times  allows  an 
aneurysm  to  begin  to  form.  The  aneurysm  is  clearly  defined  by  71  ns.  The 
velocity  oicture  at  this  time  shows  a  jet  of  high  velocity  material  con¬ 
fined  to  the  region  near  the  axis.  The  electron  number  density  contours 
are  closely  spaced  at  the  location  of  the  shell.  It  is  this  shell,  the 
locus  of  points  at  which  there  are  large  gradients  in  electron  number 
density,  that  is  to  be  compared  to  shadowqraphs  of  the  experiments.  The 
aneurvsm,  while  evident  in  Fiqure  4c,  is  not  quite  as  large  as  that  seen  in 
the  experiment  at  71  ns. 
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Figure  4a.  Velocity,  density,  temperature,  and  electron  number  density 

pictures  at  10  ns.  Flat  laser  pulse;  5  percent  preionized  1  cm 
from  the  target. 


Figure  4b.  Velocity,  density,  temperature,  and  electron  number  density 
pictures  at  30  ns.  Flat  laser  pulse;  5  percent  preionized 
from  the  target. 
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Figure  4c.  Velocity,  density,  temperature,  and  electron  number  densit, 
pictures  at  70  ns.  Flat  laser  pulse;  5  percent  preionized 
from  the  target. 


Heating  of  the  narrow  laser  channel  induces  a  weak  cylindrical 
shock  as  material  flows  radially  outward.  Consider  a  particular  location 
some  distance  from  the  target  surface.  If  the  density  at  that  axial  loca¬ 
tion  in  the  channel  falls  before  debris  from  the  laser-ablated  target 
reaches  that  location,  the  debris  will  advance  more  rapidly  in  the  axial 
direction  within  the  semi -evacuated  laser  channel.  Without  increasing  the 
total  laser  energy,  a  hotter  channel,  and  hence  a  stronger  cylindrical 
shock,  is  obtained  with  a  larger  initial  fractional  ionization  of  the  ambi¬ 
ent  medium.  Figure  5  shows  accumulated  eneroy  as  a  function  of  time  for 

the  duration  of  the  109  J  laser  pulse.  Slightly  less  than  0.1  J  from  the 

incident  laser  energy  is  absorbed  in  the  channel.  The  channel  is  here 
defined  as  that  material  with  mass  density  less  than  twice  the  ambient 
density  of  8.2  10~6  gm/cc  that  is  in  the  path  of  the  laser.  Of  course, 
simulations  with  no  preionization  (and  hence  no  absorption  of  laser  energy 
by  the  ambient  medium)  produce  spherical  blast  fronts  without  aneurysms. 

Figure  6  illustrates  electron  number  density  and  plasma  velocity 
in  the  region  close  to  the  target  on  a  finely-zoned  grid  at  2  ns  intervals 

^or  the  first  10  ns  of  the  simulation.  One  can  see  clearlv  the  formation 

of  the  laser-heated  channel  outside  the  exoandina  blast  front,  and  a  jet  of 
plasma  moving  ahead  of  the  blast  inside  the  channel. 

We  have  chosen  not  to  compute  the  effects  of  radiation  cooling  on 
the  hot  plasma  for  the  simulations  discussed  in  this  section,  so  the 
temperature  of  hot  spots  are  too  large.  This  defect  should  not  change  the 
conclusions  concerning  the  laser-heated  channel.  Our  computations  show 
that  the  laser  channel  is  not  significantly  heated  by  thermal  conduction 
from  hot  regions  near  the  target.  Thermal  conduction  from  an  essentially 
point  target  would  not  be  confined  to  a  narrow  channel,  but  would  isotro¬ 
pically  heat  the  ambient  medium. 
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Figure  5.  Laser  energy  versus  time  for  MACH2  simulation  of  NRL  85-544 
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The  issues  of  radiative  cooling  and  thermal  conduction  are 
addressed  via  numerical  simulation  in  the  next  section.  However,  we  can 
estimate  analytically  the  importance  of  thermal  conduction  over  large 
length  scales.  Conduction  of  heat  energy  from  hot  regions  obeys  a  tempera¬ 
ture  diffusion  equation 

<12) 

where  we  assume  Spitzer  conductivity,  <,  that  is  proportional  to  T5/2. 

In  addition,  thermal  flux  is  limited  to  some  fraction,  t,  times  the  energy 
flux  carried  by  free-streaming  electrons,  i.e., 

<  V  T  <  x  u  v  (13) 

where  t  is  often  taken  to  be  0.4,  u  is  the  electron  internal  energy 
density,  and  v  is  the  electron  thermal  speed.  The  Spitzer  thermal  conduc¬ 
tivity  for  a  500  eV  plasma  is  approximately  3.6  109  J/m-sec-eV.  On  length 
scales  of  1  mm,  the  time  scale  for  thermal  conduction  near  the  target  where 
the  density  is  3  10"3  gm/cc  is  60  ns.  This  is  long  compared  to  the  laser 
pulse  length,  so  thermal  conduction  in  the  hot,  moderately  dense  region 
near  the  target  is  not  expected  to  be  a  dominant  effect. 

The  debris/air  interface  outside  the  laser  channel  lies,  at  a 
particular  time,  at  the  radius  of  a  blast  wave  that  emanates  spherically 
from  a  point  source  according  to  the  similarity  relation  (Zel'dovich, 
Raizer,  1966,  Ref.  15): 


r  ( t )  =  f  (  y) 


cm 


(14) 


for  energy  (E)  in  Joules,  ambient  pressure  (P)  in  Torr,  and  time  (t)  in  ns. 
f(Y)  is  a  function  of  the  adiabatic  index  and  is  generally  of  order  0.1. 

The  NRL  experiments  fit  Equation  (14)  with  f(r)  *  0.123  (Ripin,  Stamper, 
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McLean,  1984,  Ref.  5).  For  a  109  0  laser  and  5  Torr  ambient  pressure,  the 
radius  of  the  blast  front  at  70  ns  is  1.25  cm.  To  facilitate  a  comoarison 
between  the  MACH2  simulation  and  experiment,  we  display  a  digitized  plot  of 
the  shadowgraph  for  NRL  shot  544  in  Figure  7a.  In  Figure  7b,  we  display 
the  contour  plot  of  the  electron  number  density  from  the  simulation  with 
five  percent  preionization,  drawn  to  the  same  scale  as  the  shadowgraph. 

The  MACH2-simulated  location  of  the  debris/air  interface  away  from  the  axis 
agrees  with  that  of  the  experiment.  Further,  this  location,  estimated  from 
Fiaures  7a  and  7b  to  lie  1.2  cm  from  the  target  surface,  agrees  with  that 
predicted  from  the  blast  wave  equation.  The  location  and  size  of  the 
central  aneurysm  simulated  by  MACH2  looks  remarkably  like  that  produced  by 
the  experiment.  Figure  8  is  a  superposition  of  the  digitized  version  of 
the  NRL  shadowgraph  and  the  electron  number  density  contour  plot  from 
MACH2.  We  speculate  that  numerical  diffusion  in  our  Eulerian  calculations 
together  with  the  lack  of  resolution  in  the  laser  channel  (recall  that 
there  are  but  4  cells  radially  through  the  channel)  causes  the  MACH2  simu¬ 
lated  aneurysm  to  be  flattened  at  its  tip  rather  than  to  protrude  as  far 
from  the  target  as  the  experimentally  generated  aneurysm. 


Aneurysms  were  seen  in  most  high-enerqy  exDeriments  done  with 
ambient  pressure  above  1  Torr.  Figure  9  is  a  shadowgraph  of  NRL  shot 
86-625  taken  early,  14  ns  after  peak  laser  power.  The  plastic  (CH)  target 
was  in  a  nitrogen  medium  at  5  Torr.  The  laser  pulse  was  3  ns  and 

delivered  238  J  towards  a  focal  spot  of  500  urn  radius.  The  pulse  was  ISI 
smoothed  and  produced  a  cylindrical  shock  about  the  laser  axis  as  seen  in 
the  picture.  We  simulated  this  NRL  shot  with  MACH2,  again  assuming  five 
percent  preionization  1  cm  from  the  focal  spot,  and  placed  3  computational 
cells  across  the  laser  spot,  figure  10  shows  the  electron  number  density 
from  the  simulation  at  the  same  time  as  the  shadowqraph  in  figure  9.  The 
simulation  shows  a  laser-induced  channel  similar  to  that  seen  in  the 
exneriment.  ciqure  11  is  a  superposition  of  a  digitization  of  the  shadow¬ 
graph  and  the  simulation,  drawn  to  the  same  scale. 
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RADIAL  DISTANCE  (cm)  _  RADIAL  DISTANCE  (c 


DISTANCE  FROM  TARGET  (cm) 

Igure  7a.  Digitization  of  NRL  laser/HANE  shadowgraph  for  shot  85-544  showing 
shell  and  aneurysm  at  70  and  140  ns. 


DISTANCE  FROM  TARGET  (cm) 


Figure  7b.  Contour  plot  of  electron  number  density  at  70  ns  for  MACH2  simu¬ 
lation  of  NRL  shot  85-544. 
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Figure  10.  Electron  number  density  from  MACH2  simulation  of  NRL  shot  86-625  at  14  ns. 
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Figure  11.  Superposition  of  Figures 


The  laser  axis  and  the  target  normal  coincide  in  nost  NRL  exoeri- 
ments  of  relevance  to  us.  There  are  a  few  exceotions,  however.  Figure  1? 
is  a  shadowgraph  from  NRL  shot  86-581  taken  at  59  ns  after  peak  power.  The 
beam,  carrying  95  ISI-smoothed  Joules,  strikes  the  target  30  degrees  from 
the  target  normal.  The  aneurysm  is  centered  about  the  laser  axis--not  the 
tarqet  normal. 

2.3  EFFECTS  OF  VARIOUS  PHYSICAL  PROCESSES  ON  ANEURYSM  FORMATION. 

The  purpose  of  this  section  is  to  investigate  the  effects  of  some 
of  the  physical  processes  that  occur  in  the  laser-ambient-target  system  and 
that  may  affect  the  oartially  ionized  channel  that  causes  aneurvsm  pro¬ 
duction  and  development. 

A  possible  criticism  of  our  earliest  models  of  the  laser-ambient- 
target  system  (Stel 1 ingwerf ,  Peterkin,  Longmire,  1986,  Ref.  8)  is  that  both 
the  target  and  the  ambient  consisted  of  a  single  material.  The  earliest 
models  included  a  target  and  an  ambient  medium  both  composed  of  aluminum. 
Room-temperature  gaseous  aluminum  does  not  exist  at  5  Torr,  so  our  model 
for  the  eouation  of  state  was  certainly  suspect.  Subseouent  models, 
including  those  discussed  in  Section  2.2,  included  a  nitrogen  medium  and  a 
cryroaenic  oitroqen  target.  These  models  are  reasonable  because  the 
equation  of  state  for  the  ambient  nitrogen  is  undoubtedly  correct  given  our 
assumption  of  local  thermodynamic  equilibrium,  and  the  target  material  that 
is  dynamic  in  the  simulation  is  heated  to  KeV  temperatures  early  in  the 
simulation  and  becomes  a  fully  ionized  nitroaen  plasma  and  hence  not  unlike 
aluminum  plasma.  As  target  material  expands  in  and  behind  the  blast  front, 
it  adiabatically  cools  to  tens  of  eV.  In  this  region  of  ohase  soace,  the 
differences  between  A1  and  N  may  be  noticeable.  Of  course  the  radiazion 
spectrum  is  quite  different  for  these  two  materials,  but  this  is  not  of 
conseauence  for  our  models  that  lack  soohisticated  radiation  transnort.  To 
ameliorate  the  inadequacy  of  a  nodel  with  a  sinale  material,  we  added  to 
MACH2  the  ability  to  transport  distinct  materials  across  an  Eulerian  calcu¬ 
lation  mesh.  The  differences  between  the  results  from  mul ti -mater i al  and 
single  material  simulations  are  not  significant,  so  earlier  conclusions 


28 


Figure  12.  Shadowgram  at  59  ns  of  NRL  shot  86-581.  Aneurysm  Is  aligned  with 
laser  beam  and  is  30°  from  the  target  normal. 


concerning  the  blast  front  and  the  aneurysm  remain  unchanged.  The  benefits 
of  two-material  simulations  are  two-fold.  The  approximate  position  of  the 
debris-air  interface  is  now  predicted  by  the  simulations,  and  the  inclusion 
of  distinct  materials  allows  the  attention  of  the  reader  to  be  focused  on 
the  most  important  physical  issues. 


An  additional  criticism  levied  against  our  earliest  models 
addressed  our  neglect  of  thermal  conduction.  Earlier  versions  of  MACH2 
contained  an  explicit  treatment  of  the  thermal  diffusion  process.  This 
treatment  required  the  entire  simulation  to  progress  at  the  numerical 
thermal  diffusion  timestep  if  that  timestep  was  the  smallest  in  the 
problem.  In  our  models  of  a  high-intensity  laser  interacting  with  a  solid 
target,  the  target  achieves  temperatures  in  excess  of  1  KeV.  Within  the 
Spitzer  model  for  a  fully  ionized  plasma,  the  thermal  conductivity  is  pro¬ 
portional  to  T5/2.  From  the  thermal  diffusion  Equation  (12),  we  extract 
an  explicit  timestep  for  this  process,  dttherm. 


dttherm 


p  Cv  z2  /< 


10 


-13 


sec 


(15) 


where  Cy  is  the  specific  heat  (which  we  take  to  be  approximately  0.77 
108  J/ (eV-kg ) )  when  the  temperature  is  on  the  order  of  1  KeV  (for  which 
<  -  5  1010  J/(m-sec-eV))  and  the  density  is  10_1+  gm/cc  on  a  mesh  with 
linear  cell  dimension,  z,  of  25  micron.  Computing  at  this  timestep  for  the 
10  ns  laser  pulse  would  require  105  computational  cycles.  At  a  typical 
MACH2  rate  of  150  cycles  per  CRAY  X-MP  cpu-minute,  such  a  computation  would 
require  over  10  hours  of  cpu  or  2  orders  of  magnitude  greater  than  other¬ 
wise  required  by  the  other  physical  processes  modeled  by  MACH2.  This  is 
clearly  undesirable.  Hence,  we  neglected  thermal  conduction  in  our 
earliest  models. 


To  ameliorate  this  shortcoming,  we  added  an  implicit  multi-grid 
treatment  (Brandt,  1977,  Ref.  16)  for  thermal  conduction  to  MACH2.  This 
treatment,  similar  to  our  existing  model  for  magnetic  field  diffusion. 


30 


allowed  us  to  compute  at  large  thermal  Courant  numbers.  This  implicit 
treatment  allowed  us  to  achieve  the  quasi  steady-state  solution  to  the 
thermal  diffusion  elliptic  equation  without  a  severe  computational  cost. 
The  results  of  simulations  that  include  thermal  diffusion  are  discussed 
below. 


In  addition,  we  investigated  some  other  physical  processes  that 
have  been  thought  at  various  times  to  be  of  possible  consequence  to  aneu¬ 
rysm  production.  We  have,  in  the  numerical  laboratory,  the  luxury  of 
choice  when  investigating  a  physical  process.  We  have  varied  the  degree  of 
preionization  of  the  ambient  medium,  investigated  the  effect  of  thermally 
generated  magnetic  fields,  allowed  the  low-density  ambient  medium  to 
radiate  energy  without  reabsorption,  included  a  1  Tesla  ambient  magnetic 
field,  and  varied  our  model  for  thermal  conduction.  Our  model  for  the 
thermal  conduction  includes  the  suppression  of  the  thermal  conductivity  in 
the  direction  orthogonal  to  a  strong  magnetic  field.  Thermal  energy  is 
carried  by  electrons,  and  their  mobility  in  a  direction  perpendicular  to  a 
magnetic  field  is  suppressed  as  the  electrons  gyrate  about  the  field  lines. 
The  electron  gyro-radius  is  inversely  proportional  to  the  magnetic 
induction,  so  thermal  conduction  is  most  suppressed  in  regions  of  the 
strongest  field  strength.  In  Section  5,  we  shown  that  the  magnetic  fields 
that  arise  from  thermal  energy  are  large,  on  the  order  of  tens  of  Tesla, 
near  the  edge  of  an  intense  laser  beam.  These  fields  are  anisotropic,  and 
hence,  may  effect  the  conduction  of  energy  via  electrons  from  hot  regions. 
This  anisotropy  may  effect  the  time  development  of  the  blast  front  and  the 
associated  aneurysm.  The  results  of  these  interesting  simulations  are 
discussed  below.  Discussion  of  the  effects  of  an  ambient  magnetic  field 
is  deferred  until  Section  6. 

Our  parameter  search  is  by  no  means  complete.  We  have  not,  for 
example,  investigated  with  our  two-dimensional  code  the  effect  of  ambient 
pressure  on  the  aneurysm,  nor  have  we  investigated  the  effect  of  tarqet 
thickness.  The  size  of  the  laser  focal  spot  and  the  magnitude  of  the 
incident  laser  energy  has  not  been  varied.  Our  modest  parameter  search 
includes  12  numerical  simulations  that  are  described  and  compared  below. 
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All  simulations  discussed  in  this  section  model  the  interaction 
between  a  1.06  micron  laser  beam  that  delivers  109  J  to  a  solid  A1  target 
through  a  5  Torr  ambient  medium  of  dry  air  which  is  80  percent  N2.  The 
half  convergence  angle  of  the  beam  is  2.86  degree.  The  laser  pulse  endures 
for  8  ns,  and  peaks  at  4  ns  with  an  essentially  Gaussian  temporal  profile. 
These  parameters  correspond  nominally  to  the  NRL  shot  85-544  with  the 
Pharos  II  facility,  and  are  the  same  as  those  simulations  discussed  in 
Section  2.2.  The  spatial  character  of  the  beam  is  uniform  with  a  discon¬ 
tinuous  drop  to  zero  intensity  at  the  edge  of  the  beam.  Previous  research 
(Stellingwerf,  Peterkin,  longmire,  1986,  Ref.  8)  indicated  that  a  spatially 
non-uniform  laser  beam  does  not  strongly  affect  the  character  of  the  aneu¬ 
rysm  in  the  absence  of  a  model  for  thermally  generated  magnetic  fields. 

The  attributes  of  the  twelve  simulations  discussed  in  this  section  are 
tabulated  below.  For  each,  the  computational  region  measures  1/2  cm  x 
1/4  cm. 


Table  1.  Parameters  that  denote  the  various  physical 
processes  for  the  12  models  discussed  in 
Section  2.3. 


Model  Fractional  Thermal  Radiation 

Number  Ionization  Conduction  Cooling  (1)  sgmf  (2)  Other 
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n 

2 

0.25 

% 

3 

1 
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3 
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5 
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7 

5 

% 

8 

5 

of 

•i 

9 

1 

%  (4) 

in 

5 

% 

li 

5 

x 

12 

5 

% 

on  off 

on  off 

oT>-  off 

on  off 

on  off 

on  (3)  off 

off  off 

on  on 

on  off 

on  off 

anisotropic  (5)  off 
isotropic  off 


off 

off 

off 

off 

off 

off 

off 

off 

off 

off  no 

target  (5) 
on 
on 
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Table  Notes: 


(1)  Our  tirade!  for  radiation  cooling  involves  no  radiative  trans¬ 
port,  but  rather  assumes  an  optically  transparent  medium 
whereby  energy  is  removed  from  the  plasma  in  proportion  to 
an  average  Planck  opacity  times  the  fourth  power  of  the 
plasma  temperature.  As  such,  it  represents  an  upper  bound 
to  local  energy  loss  via  radiative  recombination. 

(2)  Self  generated  magnetic  fields  (sqmf)  arise  by  conversion  of 
thermal  energy  in  the  plasma  to  magnetic  energy  in  regions 
where  the  pressure  and  density  gradients  are  not  parallel : 
i.e.,  near  the  edae  of  the  laser  beam.  Self-generated 
maqnetic  fields  are  discussed  in  detail  in  Sections  a  and  S. 

(3)  The  implicit  model  for  thermal  conduction  was  made  more 
accurate  for  this  simulation  by  limiting  the  timestep  to 
10"11  sec  per  computational  cycle.  This  is  the  only  differ¬ 
ence  between  simulations  5  and  6.  There  is  somewhat  more 
thermal  diffusion  in  model  5  than  in  the  companion  model  5. 

(4)  In  all  but  model  9  above,  our  preionization  of  the  ambient 
medium  falls  as  1/r2  from  the  laser  focal  spot  and  has  the 
value  quoted  in  the  second  column  above  at  a  point  1  cm  from 
the  laser  focus.  For  model  9,  the  preionization  still  falls 
as  1/r2,  but  from  points  1  cm  distant  from  the  center  of  the 
laser  focus  in  the  direction  tangent  to  the  target  surface. 
This  is  oerhaps  the  best  that  we  can  do  with  our  2  1/2 
dimensional  code  in  model ina  a  laser  beam  that  misses  the 
taraet,  but  interacts  with  the  blast  wave  (and  a  preionized 
medium)  that  is  produced  bv  a  second  laser  beam  that  does 
strike  a  solid  target,  traveling  parallel  to  the  pirst  at  a 
distance  1  cm  from  the  first.  This  simulation  could  be 
compared  to  the  NRL  Pharos  III  experiment  86-390. 

(5)  Model  10  includes  the  1/r2  model  of  preionization,  but 
contains  no  target--clearly  an  unphysical  situation.  How¬ 
ever,  this  model  is  of  interest  because  it  iPustrates  the 
behavior  of  the  laser-induced  cylindrical  shock  without  its 
interaction  with  a  spherical  blast  wave. 

(6)  "he  issue  of  the  isotropv  of  thermal  conduction  is  relevant 
only  in  the  presence  of  a  strong  magnetic  field  (suf¬ 
ficiently  stronq  so  the  product  of  the  self-collision  time 
and  the  gyro-frequency,  which  is  linear  in  the  magnetic 
induction,  is  greater  than  unity).  Onlv  models  U  and  12 
contain  magnetic  fields  of  any  kind,  and  these  are 
thermal! y-generated  from  thermal  enemy  in  the  plasma. 


Fiqures  13  through  13  are  snapshots  of  various  physical 
quantities  of  interest  in  the  models  listed  aoove  at  the  time  of  10  ns — 

2  ns  after  the  laser  pulse  ends.  To  facilitate  a  comoarison  between  the  12 
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gure  15. 


Material  contours  noting  the  approximate  location  of  the 
debris-air  interface  for  the  12  simulations  in  Section  2.3 


Fiqure  18.  Thermal  enerqy  flux  vectors  for  the 
Section  2.3. 


models,  each  figure  contains  12  snapshots  of  a  single  physical  attribute 
from  the  simulations — one  from  each  model.  Recall  that  the  axis  of  symme¬ 
try  is  a  vertical  line  on  the  left  of  each  snapshot.  The  laser  enters  from 
above  and  is  constrained  to  the  channel  that  is  marked  by  a  line  to  the 
immediate  right  of  the  symmetry  axis. 

Figure  13  illustrates  velocity  vectors  and  the  magnitude  of  the 
larqest  vector  is  noted  below  each  snaoshot.  As  we  increase  the  degree  of 
oreionization  in  models  1  through  5,  the  size  of  the  orotrusion  in  the 
laser  channel  qrows.  This  is  the  aneurysm,  and  the  distance  that  it 
extends  beyond  the  surrounding  spherical  blast  front  is  sensitive  to  the 
degree  of  preionization.  The  maximum  velocity  in  models  1  through  5  is  the 
same:  5.7  107  cm/sec.  Model  6  is  most  similar  to  model  5;  the  accuracy  of 
the  implicit  thermal  diffusion  model  is  somewhat  better  in  model  6.  The 
velocity  snapshots  look  similar,  but  the  maximum  velocity  is  somewhat 
smaller  in  model  6  indicating  that  the  increased  thermal  smoothing  sup¬ 
presses  only  the  maximum  plasma  speed.  Model  7  is  5  percent  preionized, 
but  has  no  thermal  conduction.  The  peak  velocity  and  the  extent  of  the 
aneurysm  are  not  as  great  as  any  of  the  models  that  are  also  5  percent  pre¬ 
ionized  but  include  thermal  conduction.  Model  3  radiates,  but  the  effect 
of  radiative  cooling  on  the  hydrodynamics  seems  from  Figure  13  not  "p  be 
great.  Model  9  shows  no  aneurysm  when  the  source  of  preionization  is  1  cm 
off-axis.  While  not  conclusive,  this  does  indicate  that  experimentally  one 
should  not  exoect  to  see  a  laser-al iqned  aneurysm  when  the  laser  beam  ioes 
not  strike  a  target  (as  is  the  case  in  NRL  shot  36-590).  Model  10  contains 
no  target,  but  incudes  the  1/r2  model  for  oreionization  that  oriainates 
from  the  laser  focal  soot  (an  unohysical  situation).  The  oeak  velocity 
profile  indicates  expansion  from  the  heated  channel.  Models  11  and  12 
include  simulations  of  thermally-generated  magnetic  fields.  The  peak 
velocities  in  these  two  models  are  suppressed  somewhat  by  the  maonetic 
fields,  but  the  relative  position  of  the  spherical  blast  front  and  tne 
ineurvsm  are  similar  to  model  5  that  has  the  same  value  for  the  oreioni¬ 
zation. 
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Figure  14  illustrates  iso-temperature  contours  on  a  logarithmic 
scale.  The  location  of  the  maximum  temperature  is  indicated  by  a  olus  (+) 
sign,  and  this  peak  value  is  noted  at  the  top  of  each  frame.  An  average 
temperature  that  is  characteristic  of  the  shell  is  noted  also  in  each 
frame.  The  temperature  contours  in  model  1  show  a  laser-heated  cnannel  not 
to  be  present,  of  course,  as  there  is  no  preionization  to  seed  laser 
absorption.  A  small  dimple  in  the  blast  front  appears  in  the  laser  channel 
at  the  smallest  nonzero  value  of  preionization  in  model  2.  As  the 
fractional  oreioni zation  increases  between  models  2  and  5,  the  cylindrical 
laser-heated  channel  becomes  more  pronounced.  Recall  chat  it  this  laser- 
induced  cylindrical  structure  that  we  interoret  as  the  oroqenitor  of  the 
aneurysm.  Interestinql v,  models  1  through  5  have  the  same  model  for 
thermal  conduction,  and  all  5  models  produce  the  same  shell  temperature  and 
maximum  temperature,  independent  of  the  fractional  oreioni  zation .  Modal  o 
is  somewhat  cooler  than  the  first  5  models  because  that  model  encourages 
thermal  diffusion.  Model  7  has  no  thermal  conduction  and  is  about  20  per¬ 
cent  hotter  than  models  1  through  5.  However,  the  relative  size  of  the 
cylindrical  laser-heated  channel  and  the  spherical  blast  front  is  the  same 
in  models  5,  6,  and  7.  This  strengthens  our  claim  that  thermal  liffusion 
is  not  an  important  mechanism  for  aneurysm  growth  at  this  moderate  value  of 
the  laser  energy.  Model  3  has  the  same  node!  for  thermal  conduction  as  the 
first  5,  but  has  our  simple  radiation  cooling  mechanism  turned  on.  'he 
cylindrical  laser-heated  channel  persists,  and  the  plasma  temperature  is 
about  10  oercent  smaller  than  the  similar  oon->*adiative  models.  Because 
the  source  of  preionization  is  further  away  from  the  laser  channel  in  model 
9,  the  laser-heated  channel  is  not  as  pronounced  as  in  model  3.  Model  10, 
without  a  target  but  with  the  same  model  of  preionization  and  thermal  con¬ 
duction  as  model  5,  shows  a  strong  laser-induced  cylindrical  structure. 
Models  11  and  12  are  5  percent  preionized  just  as  model  5,  but  include 
thermally-qenerated  magnetic  fields,  'he  structure  of  the  ’aser-heated 
channel  is  similar  for  these  3  models.  Notice,  however,  that  nodel  11, 
with  anisotropic  thermal  conduction,  produces  Varner  temperatures  in  the 
shell  and  in  the  bubble. 
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Figure’ 15  illustrates  electron  number  density  iso-contours  on  a 
loqarithmic  scale.  The  structure  of  these  snapshots  essentially  mimics  the 
corresponding  temperature  features  of  Figure  14. 

Figure  15  illustrates  the  approximate  position  of  the  deDris-air 
interface.  The  broadness  of  the  band  is  a  consequence  of  numerical  dif¬ 
fusion  of  material  across  our  Eulerian  mesh.  The  middle  of  this  band 
represents  the  location  of  the  interface.  This  interface  is  essentially 
spherical  outside  the  laser  channel.  Models  1  throuah  5  show  that  the 
debris-air  interface  protrudes  ahead  of  the  blast  front  at  Doints  within 
the  channel  by  an  amount  proportional  to  the  fractional  preionization. 

This  indicates  that  the  aneurysm  is  made  of  target  material.  The  laser- 
heated  channel  exists  far  in  front  of  the  aneurysm  and  consists  of  ambient 
material.  Model  ID  has  no  target,  hence  no  debris. 

Iso-density  contours  on  a  logarithmic  scale  are  illustrated  for 
our  12  models  in  Figure  17.  They  clearly  illustrate  the  position  of  the 
spherical  blast  front  and  the  1 aser-induced  cylindrical  blast  wave.  ”he 
nass  density  within  the  channel  is  less  than  the  ambient  channel,  and  :hi$ 
allows  target  debris  to  exoand  preferenti ally  along  the  laser  axis  and 
create  the  aneurysm. 

Finally,  Figure  13  illustrates  thermal  flux  vectors  for  10  of  the 
12  models.  Recall  that  model  7  does  not  compute  thermal  conduction.  In 
all  except  model  11,  the  thermal  flux  is  directed  away  from  the  hot  shell 
in  an  essentially  symmetric  manner.  However,  in  model  11,  the  thermal  fux 
is  seen  mostly  within  the  laser  channel,  and  is  suppressed  outside  the 
channel  where  thermally-generated  maqnetic  fields  are  strong.  These  fields 
exist  mostly  within  the  shell  away  from  the  symmetry  axis,  and  are  toroidal 
in  character.  (Details  are  qiven  in  Section  5.)  It  is  possible  that  for 
low  energy  laser  pulses,  where  the  preionization  is  not  as  great  as  we  have 
assumed  in  these  models,  the  anisotropic  thermal  conduction  mav  aniso- 
tropize  the  blast  wave. 
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2.4 


CONCLUSIONS. 


Simple  analytic  analysis  suggests  that  the  channel  of  ambient 
room  temperature  air  through  which  the  laser  propagates  can  be  heated  to 
many  eV  via  inverse  bramsstrahlung  laser  absorption  by  that  air  if  it  has  a 
non-zero  ionization  fraction.  The  chanqe  in  temperature  of  the  ambient 
gas  is  proportional  to  the  ionization  fraction  of  the  gas  (Equation  7). 
Therefore,  there  must  exist  some  mechanism  for  preionizing  the  room  temper¬ 
ature  ambient  medium  in  order  to  initiate  laser  absorption.  There  is  some 
experimental  evidence  that  the  ambient  medium  is  oreionized  by  a  few  oer- 
cent  at  a  distance  of  1  cm  from  the  target  surface,  perhaps  due  to  X-rays 
emanatinq  from  the  1 aser-irradiated  target  and  other  hot  plasma.  This 
degree  of  preionization  is  sufficient  to  initiate  heating  of  the  channel 
material  by  the  laser. 

MACH2  simulations  give  conclusive  evidence  that  a  laser-heated 
channel  leads  to  the  formation  of  a  central  aneurysm  in  NRL  laser/HANE 
experiments.  The  seed  for  the  formation  of  an  aneurysm  aooears  by  the  time 
if  peak  laser  oower  as  a  weak  cylindrical  shock  centered  aoout  the  laser 
axis.  If  the  heating  is  appreciable  far  from  the  target  on  a  time  scale 
:omparable  to  the  lenath  of  the  laser  oulse,  an  aneurysm  will  oe  evident  at 
’ater  times.  Heating  ’n  the  'aser  channel  is  proportional  to  the  ioni¬ 
zation  fraction,  as  .veil  as  to  the  laser  intensity,  so  the  aneurysm  nay 
scale  with  laser  energy  (this  is  consistent  with  snadowqraDhs  from  low 
energy  shots).  Figure  3,  a  superposition  of  the  electron  number  density 
generated  by  MACH?  for  a  simulation  that  is  five  oercent  oreionized  :m 
*rom  the  target  surface  upon  a  digitization  of  an  URL  shadowgraph  from 
experimental  shot  number  35-544  illustrates  excellent  agreement  between 
experiment  and  theory.  cigure  11,  a  superposition  of  an  early-time  shadow¬ 
graph  from  URL  exoeriment  36-525  and  a  picture  from  a  computer  simulation 
of  that  shot  at  the  same  time,  shows  the  laser-induced  cylindrica1  shock 
front  from  the  exoeriment  and  a  comouter  simulation  of  the  same. 
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SECTION  3 


A  SIMPLE  MOOEL  FOR  LASER  CHANNEL-HEATING  IN  FEW-TORR  PLASMA 

A  hot,  low  density  laser  channel  has  been  shown  to  be  imoortant 
in  many  of  the  simulations  discussed  in  this  report.  The  cylindrical  blast 
*rom  the  heating  is  seen  in  the  shadowgraphs  at  early-time  (10  ns),  and  we 
has  shown  that  it  is  the  primary  cause  of  the  on-axis  aneurysm  seen  in  the 
NRL  experiments.  Because  of  this  direct  cause-and-effect  rel ationship,  it 
is  of  interest  to  inquire  how  the  laser  heating  depends  on  the  experimental 
parameters.  Since  this  heating  takes  place  during  the  laser  pulse  (1  - 
5  nsl,  it  is  reasonable  to  ignore  the  hydrodynamics  and  other  slower 
orocesses,  and  consider  the  extent  and  nature  of  the  initial  oeated  region 
as  a  separate  oroblem.  In  this  section  we  develop  a  simple  model  of  the 
heated  channel  that  aarees  very  well  with  the  detailed  two-dimensional 
models,  and  allows  computation  of  scaling  laws  for  the  heatinq  process. 

'•/e  will  thus  demonstrate  that  direct  inverse-bremsstrahl ung  deposition  of 
laser  energy  can  account  for  the  ooserved  and  computed  features  of  the 
channel,  and  that  no  shock  or  detonation  mechanisms  appear  to  be  ooerafnn 
in  this  regime  durina  the  laser  oulse.  "he  model  used  here  is  a  direct 
extension  of  that  of  Section  3.1,  used  to  motivate  the  channel -heating 
simul ations. 

"he  incident  ’aser  enemy  orooaaates  through  the  underdense 
ambient  medium  and  deposits  the  bulk  of  its  energy  at  the  critical  surface, 
whic*-  is  the  target  surface  at  the  onset  of  the  laser  pulse.  Exper ■’'mental 
evidence  indicates  that  a  orief  ourst  of  X  rays  are  emitted  from  the  target 
during  the  early  stages  of  the  laser  pulse,  "hese  X  rays  ohotoionize  the 
ambient  gas  (bound-free  absorption)  and  orovide  the  seeds  for  inverse 
bremsstrahl unqen  (free-free)  absorotion  of  the  incident  laser  beam  bv  the 
ambient  gas.  We  do  not  attempt  to  model  this  initial  target  heatinq  and 
ohotoioni zat'on  process  here,  out  instead  simol  •/  assume  that  the  ambient 
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gas  is  initially  ionized  with  an  ionization  fraction  that  drops  as  1/r2 
where  r  is  the  distance  from  the  target.  The  degree  of  ionization  is  found 
experimentally  to  be  in  the  range  of  one--few  percent  at  a  distance  of  1  cm 
from  the  target. 

The  laser  intensity  in  the  channel  is  given  by  I  =  I0  exp(-/  K 
dr),  or,  for  small  absorption,  I  =  Io  (1-/  K  dr),  where  I0  is  the  incident 
intensity,  and  K  is  the  absorption  coefficient  (units  of  nr1).  The  formula 
for  the  inverse  bremsstrahl ung  absorption  coefficient  has  been  given  by 
Johnston  and  Dawson,  1973  (Ref.  17).  For  the  present  case  (electron 
density  much  less  than  critical),  this  coefficient  is  given  by  Equation  (3) 
and  is  repeated  for  clarity: 


<  * 
V 


a  b  A 


2  2, 

o  i  n  a 

Tzm — 


(16) 


where  X  is  the  laser  wavelenqth,  ln(A)  is  the  coulomb  logarithm  (varies 
from  1-5  for  these  cases),  '  is  the  number  of  free  electrons  per  nucleon, 
=  maxU,  1),  4  is  the  atomic  number  of  the  gas,  p  the  density  of  the 
qas,  and  T  the  local  gas  temperature.  The  product  3b  is  1.5666  *  1018/m 
for  mks  units  with  T  in  eV.  The  important  point  here  is  that  the  absorp¬ 
tion  depends  strongly  on  ionization  state  and  inversely  on  the  ambient 
temperature.  The  factor  accounts  for  the  presence  of  neutral  ambient 
atoms  at  low  degrees  of  ionization. 

The  time  derivative  of  temperature  at  a  point  in  the  channel  is 

given  by: 


dT  _  %  I 
"3TT  ~  p  uv 


(17) 
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where  <  is  given  by  Equation  (16),  I  is  the  local  laser  intensity,  and  1 

v  v 

is  the  specific  heat.  The  intensity  will  depend  on  the  time  variation  of 
the  laser  energy  as  well  as  the  distance  from  the  tarqet.  For  a  oerfect 
gas,  Cv  =  1.44  *  108  (1  +  c)/A  J/eV-kg. 

We  may  integrate  Equation  (17)  directly  to  obtain  an  estimate  of 
the  final  temperature  of  the  laser-heated  channel.  We  assume  fully  ionized 
nitrogen  as  the  ambient  gas,  p  =  5  Torr,  A  =  1  micron,  1 n ( A )  =  3,  and  find 

Tf  -  (6  E/r2)2/5  ,  (13) 

where  E  is  the  laser  energy  in  Joules,  and  r  is  the  local  channel  radius  in 
centimeters,  T,  the  final  temperature  in  eV,  and  we  have  assumed  that  the 
initial  channel  temperature  is  small.  Given  these  assumptions.  Equation 
(18)  in  consistant  with  Equation  (7).  Taking  E1  =  100  J,  and  r  =  0.05  cm, 
we  find  Tf  =  142  eV,  so  a  hot  channel  is  expected  and  should  extend  a  con¬ 
siderable  distance  alonq  the  laser  channel,  since  the  dependence  on 
distance  is  weak.  On  the  other  hand,  the  final  temoerature  is  aoDroxi- 
mately  1  inearl v  dependent  uoon  the  deqree  of  ionization,  oroducina  a  feed¬ 
back  mechanism  that  should  shut  off  the  heating  at  a  distance  where  the 

available  deoosited  energy  (s  too  small  to  ionize  the  ambient  oas. 

1  better  estimate  of  the  channel  temoerature  as  a  function  of 
time  (during  the  laser  pulse)  and  distance  from  the  target  is  obtained  by 
numerical  integration  of  Equation  (17)  with  time  dependent  quantities 
■ncluded.  cor  this  estimate  a  Gaussian  pulse  shape  is  used  tor  the  1 aser 

lower,  oeakina  at  i  ns,  endinn  at  3  os,  the  nosition  in  the  channel 

determines  the  level  of  preionization  and  the  channel  radius,  the  Goulomb 
logarithm  is  given  by  (Spitzer,  1962,  Ref.  13): 

ln(A)  =  'n (min  (4.9e9  ~L’;>  A/c,  8.08e8  T  xi) 
with  Inf  A)  >=  1. 
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The  ionization-  fraction,  5,  is  given  by  a  fit  to  Saha  equilibrium 
at  5  Torr  pressure  of  atomic  nitrogen  shown  in  Figure  19.  At  this  pressure 
nitrogen  begins  to  ionize  at  about  1  eV,  is  five  times  ionized  at  10  - 
30  eV,  and  is  fully  ionized  at  90  eV. 

The  specific  heat,  C  ,  for  nitrogen  is  shown  in  Figure  20.  The 
irregular  line  is  taken  from  a  tabular  equation  of  state  for  nitrogen  at  5 
Torr  pressure,  the  line  labeled  "fit"  is  the  rough  fit  used  for  the  present 
models,  and  the  dashed  line  is  the  oerfect  gas  value  ignoring  ionization 
energy.  By  comparing  the  areas  above  and  below  the  dashed  line,  we  con¬ 
clude  that  for  temperatures  in  the  0.1  -  100  eV  range  most  of  the  laser 
energy  absorbed  in  the  channel  is  needed  to  ionize  the  gas. 

We  have  integrated  Equation  (17)  for  target  distances  ranging 
from  0.1  to  3  cm.  We  call  the  result  of  these  integrations  the  "zero 
dimensional"  (or  0  0)  model,  since  each  point  in  the  channel  is  heated 
completely  independently  of  all  other  points.  The  results  depend  on  the 
degree  of  preionization  of  the  channel.  Figure  21  shows  the  case  of  5  per¬ 
cent  ionized  at  1  cm  from  the  target.  Temperatures  are  shown  as  a  function 
of  distance  along  the  channel  for  times  of  0  (room  temperature)  to  8  ns 
fend  of  laser  pulse).  The  result  from  Equation  (18)  is  labeled  "simple 
estimate."  As  expected.  Equation  (18)  is  a  good  estimate  of  channel  heat¬ 
ing  for  regions  that  become  fully  ionized  (0  -  0.5  cm),  but  beyond  0.5  cm 
the  temperature  in  the  channel  drops,  with  little  heating  beyond  0.9  cm. 

An  ionization  front  moves  to  larger  radii  in  this  figure  with  a  mean 
velocity  of  about  108  cm/s,  heating  the  gas  from  ambient  temperatures  to 
about  100  eV,  but  it  is  not  a  progressive  wave,  since  no  coupling  between 
different  radii  is  included  in  the  model. 

The  result  is  a  channel  heated  to  about  100  eV  extending  from  the 
target  to  about  0.9  cm  along  the  channel.  Tip's  heated  region  expands 
hydrodynamical ly  as  a  cylindrical  shock  wave  and  forms  a  low  density 
channel  at  later  times.  The  spherical  blast  wave  expands  more  rapidly  in 
this  channel,  resulting  in  aneurism  formation.  When  the  blast  wave  moves 
beyond  the  heated  region,  the  aneurism  should  decay. 
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Figure  19.  Fit  to  tabular  equation  of  state  for  the  ionization  fraction  of 
nitrogen  at  5  Torr  pressure  as  a  function  of  temperature.  These 
numbers  agree  with  Saha  equilibrium  values. 
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Figure  20.  Fit  to  tabular  equation  of  state  for  the  specific  heat  of  nitrogen 
at  5  Torr  pressure  as  a  function  of  temperature.  Both  the  exact 
tabular  values  (irregular  solid  line)  and  the  simple  approximation 
used  here  are  shown.  The  dashed  line  is  an  estimate  of  the  perfect 
gas  value  obtained  when  ionization  energy  is  ignored. 
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Figure  21.  Run  of  temperature  versus  position  in  the  channel  at  various  times 
during  the  laser  pulse  for  case  1  of  Table  2.  Also  shown  is  the 
result  obtained  from  Equation  (18)  labeled  "simple  estimate". 


Figure  22.  Temperature  at  the  head  of  the  heated  channel  as  a  function  of  the 
convergence  angle  of  the  laser  beam.  The  two  laser  configurations 
used  at  NRl  are  shown. 
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Table  2  shows  the  results  of  a  number  of  models  run  with  various 
values  for  the  laser  energy  (E),  ambient  pressure  (P),  and  degree  of  pre¬ 
ionization  at  1  cm  ( Co ) •  The  listed  parameters  are:  L  =  extent  of  heated 
region  from  target  of  preionization  of  the  ambient  gas  (cm),  T^  ^  = 
temperature  at  the  top  of  the  ionization  wave  (eV),  T.|  =  temperature  at 

the  bottom  of  the  ionization  wave  (distance  =  size)(eV). 


Table  2.  Laser  heated  channel  characteristics. 


E 

X 

P 

'*0 

angle 

L 

Thigh 

^1  ow 

1 

109 

1.0 

5.0 

5.0 

2.86 

0.9 

25 

1.0 

2 

109 

1.0 

5.0 

2.5 

2.86 

0.7 

45 

1.0 

3 

109 

1.0 

5.0 

10.0 

2.86 

1.1 

20 

1.0 

4 

109 

1.0 

2.5 

5.0 

2.86 

0.8 

16 

1.0 

5 

109 

1.0 

10.0 

5.0 

2.86 

1.0 

40 

1.0 

6 

50 

1.0 

5.0 

5.0 

2.86 

0.8 

16 

1.0 

7 

200 

1.0 

5.0 

5.0 

2.86 

1.0 

30 

1.0 

8 

109 

1.0 

5.0 

5.0 

10.0 

0.6 

45 

1.0 

9 

109 

1.0 

5.0 

5.0 

8.0 

0.65 

22 

1.0 

10 

109 

1.0 

5.0 

5.0 

6.0 

0.7 

20 

1.0 

11 

109 

1.0 

5.0 

5.0 

1.5 

1.1 

50 

1.0 

12 

109 

1.0 

5.0 

5.0 

0.0 

2.5 

400 

1.0 

13 

109 

0.5 

5.0 

5.0 

2.36 

0.7 

20 

1.0 

14 

109 

2.0 

5.0 

5.9 

2.36 

1.2 

100 

1.0 

The 

size 

of  the 

heated 

region  in 

models 

1  -  7 

is  fit 

to  an 

accuracy  of  2  percent  by  the  scaling  law: 


L  =  n. 19  En,Ifi  P0,16  cn0,33 


(20) 


where  L  is  the  length  of  the  heated  region  in  centimeters,  P  is  the 
ambient  pressure  of  N2,  and  c0  is  the  percent  preionization  at  1  cm 
crom  the  target.  For  example,  for  c0  =  5%,  E  =  100  .),  and  0  =  5  Torr, 
we  find  L  =  0.80  cm. 
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The  scaling  law  for  a  simple  blast  wave  (Equation  (141)  is: 

R  =  0.123  (E/P)0-2  t0*4  .  (21) 

The  energy  scaling  of  Equation  (20)  is  that  of  a  normal  blast  wave,  but  the 
pressure  scaling  is  the  opposite  of  the  blast  wave  case.  This  means  that, 
as  the  ambient  pressure  is  increased,  the  size  of  the  aneurism  produced 
increases  'while  the  position  of  the  blast  wave  at  a  fixed  time  decreases, 
causing  the  pressure  deoendence  to  be  the  most  obvious  in  the  experiments. 
Increasing  the  laser  energy  causes  both  the  heated  channel  and  the  blast 
wave  position  to  increase  similarly,  so  the  size  of  the  aneurism  appears  to 
be  independent  of  this  quantity. 

The  results  from  Table  2  for  the  temperature  at  the  head  of  the 
ionization  wave  are  fit  to  10  percent  accuracy  by  the  scaling  law: 

t  _  o  a~  j0.48  „0.o6  -0.53  ?  ?? ) 

Thigh  -  2-4°  ?  "0  • 

""he  temperature  increases  with  laser  energy  and  ambient  pressure,  out 
decreases  strongly  as  the  oreionization  is  increased.  The  effect  of  /dry¬ 
ing  the  preioni zation  is  simply  to  change  L  while  maintaining  the  tempera¬ 
ture  as  a  function  of  distance  nearly  unchanged,  thus  producing  higher 
"head"  temperatures  for  the  shorter  channels. 

The  models  3-12  in  Table  2  show  the  dependence  uoon  the  con¬ 
vergence  angle  of  the  laser  beam.  We  see  that  for  longer  focal  length 
lasers  the  size  of  the  aneurysm  is  ouch  longer,  while  the  temoerature 
increases  for  both  small  and  large  beam  angles.  This  probably  explains  why 
the  aneurysm  became  much  more  obvious  when  the  laser  at  NRL  was  upgraded  to 
a  higher  enerqv  configuration  and  the  focal  length  was  changed  from  a  P/5 
'4.75  degree  half-angle)  to  a  f/10  (2.35  degree  half-anole)  focal  ratio. 

The  dependence  of  size  (L)  on  half-anole  of  the  beam  ':s  given  approximate  1  y 
by 
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L  =  1.26  (half-angle) 


-0.32 

where  the  angle  is  in  degrees  and  L  is 
ture  upon  angle  is  shown  in  Figure  22. 
the  zero  angle  model  has  been  omitted, 
case. 


,  (23) 

in  cm.  The  dependence  of  tampera- 
In  both  Equation  (23)  and  Figure  22 
and  should  be  regarded  as  a  limiting 


We  now  consider  the  details  of  the  absorption  process.  Figure  23 
shows  the  absorption  coefficient,  ie  (Equation  (16))  versus  distance  for 
model  1  in  Table  2  at  times  1  -  3  ns.  At  very  early  times  the  absorption 
occurs  near  the  target,  but  as  early  as  2  ns  the  peak  absorption  has  moved 
away  from  the  target  and  henceforth  occurs  near  the  head  of  the  channel 
where  the  number  of  free  electrons  is  high  and  the  temperature  is  still 
relatively  low.  Absorption  near  the  target  is  very  low  at  later  times  due 
to  the  very  high  temperatures  in  this  region.  Integrating  these  curves,  we 
obtain  the  total  fraction  absorption  versus  time  curve  shown  in  Figure  ?A. 
The  absorotion  increases  during  the  first  half  of  the  laser  pulse  (0  - 
a  ns),  then  levels  out  at  about  0.5  percent.  This  results  in  about  0.5  J 
energy  absorbed  in  the  channel  for  a  100  J  laser  pulse,  in  agreement  with 
the  two-dimensional  models. 

Finally,  models  13  and  14  show  the  effect  of  varying  the  laser 
wavelenqth.  The  size  of  the  heated  region  is  given  to  an  accuracy  of  1 
percent  by 

n  qq 

L  =  0.91  X  ,  (24) 

where  L  is  in  cm  and  x  is  in  microns.  'rhe  temperature  of  the  plume  is 
•'elatively  constant  below  x  =  1  micron,  and  increases  rapidly  for  longer 
wavelengths  to  100  eV  at  2  microns  and  several  keV  at  10  microns  (not 
shown) . 


This  model  makes  no  allowance  for  a  variety  of  possible  coolina 
mechanisms,  and  thus  represents  an  upper  limit  estimate  of  the  channel 
heating.  In  the  first  3  ns  of  target  irradiation,  however,  there  is  no 
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ABSORPTION  (I) 


.2 


0.2  0.4  0.6  0.8  1.0  1.2 

DISTANCE  FROM  TARGET  (cm) 

Figure  23.  Absorption  coefficient  (Equation  (16))  versus  position  in  the 
channel  at  various  times  during  the  laser  pulse  for  model  1  of 
Table  2. 


TIME  (ns) 


Figure  24.  Instantaneous  fractional  absorption  as  a  function  of  time  during 
the  laser  pulse  for  model  1  of  Table  2.  Peak  laser  power  occurs 
at  4  ns. 


time  for  hydrodynamic  expansion,  thermal  conduction,  or  radiation  cooling, 
so  these  estimates  should  be  quite  reasonable,  ^n  fact,  the  two- 
dimensional  simulations  show  almost  identical  temperature  profiles  to  those 
shown  here,  proving  that  this  model  is  an  accurate  representation  of  the 
laser  channel  heating. 
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SECTION  4 


CAN  THERMALLY  GENERATED  MAGNETIC  FIELDS 
8E  INCORPORATED  INTO  A  SINGLE  FLUID  MODEL? 

4.1  INTRODUCTION. 

This  section  addresses  the  question  of  whether  it  is  reasonaole 
to  include  the  process  known  as  thermally  (or  self)  generated  magnetic 
fields  into  a  single  fluid  magnetonydrodynamics  (MHD)  model.  Our  con¬ 
clusion  is  that  it  is  indeed  possible  to  model  this  effect  in  this  way.  An 
outline  of  the  path  along  which  we  travel  to  arrive  at  this  conclusion 
follows.  We  then,  in  the  following  section,  illustrate  the  thermal  gener¬ 
ation  of  magnetic  fields  with  a  two-dimensional  MHD  code,  MACH2. 

4.2  THEORY  FOR  EACH  PLASMA  SPECIES. 

We  begin  with  the  fluid  equations  for  charged  particles  of  spec¬ 
ies  s.  The  species  are  either  ions  (s  =  i)  or  electrons  fs  =  e),  each  of 
which  is  characterized  by  a  mass,  n  and  charge,  a.,  "he  momenta  of  the 
particles  are  assumed  to  be  influenced  only  by  electric  and  magnetic 

fields,  E  and  R,  and  by  collisions,  "he  eouation  of  motion  for  an  an sample 

■> 

of  particles  of  species,  s,  with  a  number  density,  n  (r,t),  follows  from 
the  Vlasov  equation  if  we  assume  a  suitably  smooth  thermal  spread  in  the 
velocity  distribution: 


m 


sns 


dvs 


Vs 


+  7P$  -  K$  =  0 


(25) 


We  will  employ  cgs  units  throughout  this  section.  The  fluid  velocity,  v  , 
is  an  ensemble  average  of  single  particle  velocities.  The  pressure,  P  ,  is 
related  to  the  temperature  by  an  equation  of  state.  The  temperature  is 
simply  a  measure  of  the  mean  sauare  thermal  speed,  V,  for  particles  of 
mass,  m  .  The  final  term  in  Equation  (25)  is  a  source  of  momentum  per  unit 
volume  that  is  caused  by  inter-soecies  collisions.  Let  us  assume  K3=  -K, 
for  a  nonviscous  fluid. 
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The  fluid  of  each  species  obeys  a  continuity  equation 


It  ("s>  *  Vvi)  ■ 0  •  <26> 

if  the  total  number  of  particles  of  each  species  is  conserved.  We  employ 
the  summation  convention  whereby  a  sum  from  1  to  3  over  repeated  indices  is 
impl ied. 

4.3  MHD  EQUATIONS. 

A  single  fluid  comprised  of  electrons  and  ions  may  be  character¬ 
ized  by  its  total  mass  density,  momentum  density,  current  density,  charge 
density,  and  thermal  energy  density  (pressure).  Each  is  defined  as  the  sum 
of  the  contributions  from  each  species.  Therefore,  with  the  convention 
e  =  4.803  10-10  esu,  the  charge  of  each  free  electron  is  qg  =  -e  and  the 
charge  on  the  average  ion  is  q  =  +?e.  ?  is  between  0  and  the  atomic  number 
of  the  ion  species.  Because  we  allow  c  less  than  unity,  we  may  include 
neutral  atoms  among  our  ions.  Hence: 


mass  density: 

p(r)  =  mini(r)  +  ma"e(r) 

(27) 

charge  density: 

°c(r)  =  e(cn.(r)  -  np(r)) 

(28) 

momentum  density: 

pv  =  mini(r)vi(r)  +  Nene( r) v0( r) 

(29) 

current  density: 

3  =  e(cni(r)vi(r)  -  ng(r)ve(r)) 

(30) 

Pressure: 

a  =  p.  +  pe  . 

(31) 
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3y  combining  the  continuity  equations  and  the  equation  of  motion 
for  each  species  with  the  above  definitions,  one  may  derive  single-fluid 
equations.  For  example,  by  summing  over  both  species  the  product  of  the 
continuity  equation  and  the  mass  for  each  species,  one  obtains  the  mass 
conservation  equation: 


ap 

3t 


+ 


V. 

J 


(pvj)  =  0 


(32) 


To  obtain  a  single-fluid  equation  of  motion,  it  is  desirable  to 
express  the  velocity  for  each  species  in  terms  of  the  center  of  mass  veloc- 
ity,  v,  and  the  net  current,  J.  One  obtains: 


(33) 


(34) 


In  order  to  proceed  further,  we  assume  quasi-neutrality:  »  0.  In  this 

case,  the  electron  number  density,  ng,  is  related  to  the  ion  number  density 
via  Equation  (28);  i.e.,  ng  =  cn^.  Hence 


and 


vi 


m 


.1  - 
e  p 


+ 

*  v 


2^1 

e  pc 


(35) 


(36) 
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Note  that  the  difference  between  the  velocities  of  the  two  species  is 
parallel  to  the  current  density. 


*  *  j 

v  .  -  v  *  — 
i  e  en 


and  the  sum  is, 


3m. 


v.  +  v  =»  2v - -  (1 

i  e  epc 


S"V 


(37) 


(38) 


neglecting  c  mg/m.  compared  to  unity,  and  using  Equation  (36),  we  find  that 
-*•  -*■ 

V  .  *  V  . 

In  going  from  the  particle  description  to  a  fluid  model,  we 
replace  the  eight  particle  variables  v. ,  v  ,  n^ ,  ng  with  the  eight  fluid 
variables,  v,  3,  p,  p  .  In  MHD,  we  restrict  our  attention  to  quasi -neutral 
fluids  with  pg  =  0.  An  additional  assumption  is  that  of  ion-electron  equi¬ 
librium,  in  which  case  the  electron  temperature  is  equal  to  the  ion  temper¬ 
ature,  T  =  T..  In  this  case,  P  =  qP  ■ ,  and  as  the  total  pressure  is 
simply  Pg  +  ?. ,  we  have  Pg  =  U/l+c)  P  and  =  (1/1+c)  P. 


By  adding  and  subtracting  the  two  particle  equations  of  motion 
Equation  (25),  we  obtain  a  fluid  equation  of  motion  and  a  generalized  Ohm's 

1  dW. 


Under  the  assumptions  of  quasi-neutral ity  and  electron-ion 
equilibrium,  we  obtain 

<39> 

as  the  eauation  of  motion  for  a  quasi -neutral  conducting  fluid.  If  fluid 
viscosity  is  included  in  our  model,  we  must  add  to  Equation  (39)  the  term 
7  •  j  where  a  is  a  viscous  dyad.  Under  these  same  assumptions,  neglecting 
terms  of  order  me/m.j  when  compared  to  unity,  we  have 
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3T 


(‘ 


m  n  c 
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X  S 


+  (l+c)meni 


7P  +  K 
m  n 
e  e 


(40) 


to  lowest  order  as  a  generalized  Ohm's  law.  The  various  current  generating 
terms  on  the  right  hand  side  of  Equation  (40)  are  the  electric  field,  E, 
the  electromotive  induction,  (v/c  x  B),  the  Hall  effect,  (J  x  B),  the 
electron  pressure,  VP,  and  a  collisional  source  term,  K.  This  source  term 
will  be  considered  to  arise  solely  from  interspecies  collisions,  and  hence 
can  be  identified  with  the  familiar  resistive  term 


en„ 

£»-—£?  (41) 

3 

where  a  is  the  electrical  conductivity  of  the  fluid.  Other  sources  of 
momentum  for  the  particles  such  as  radiation  pressure  or  a  thermoelectric 
force  are  not  considered  here.  We  may  use  the  mass  conservation  Equation 
(32)  t o  express  the  left  hand  side  of  Equation  (40),  the  electron  inertia 
term,  as 


My  -  k  $  +3H  •  m 

The  electric  and  magnetic  fields,  and  the  current  density,  are 
related  via  Maxwell's  equations.  By  solving  Equation  (40)  for  the  electric 
field  and  taking  its  curl,  we  obtain  from  Faraday’s  law  an  equation  for  the 
time  derivative  of  the  magnetic  field.  We  express  the  current  density  in 
terms  of  the  curl  of  the  magnetic  field  (Ampere's  law)  by  neglecting  the 
displacement  current.  Such  neglect  is  consistent  with  our  assumption  of 
quasi-neutrality  whereby  the  charge  density  is  approximately  zero,  as  is 
it's  partial  derivative  with  respect  to  time.  After  performing  some 
algebra,  we  obtain  the  equation  for  magnetic  field  generation: 
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(43) 


m  c 
e 

"ee 


d 

1? 


(7  x  3)  -  (7  x  J)  (7*V) 


The  current  density  occurring  in  the  last  term  can  be  replaced  with 

+  ♦  , 
c/4ir  7  x  b,  and  the  curl  of  J  with  -c/4 w  7  B. 

Note  that  if  the  gradient  of  the  number  density,  the  velocity, 
and  the  gradient  of  the  conductivity  all  vanish,  and  if  we  can  ignore  the 
d(v  x  j)/dt  term,  we  recover  the  familiar  diffusion  eauation  for  the 
magnetic  field: 


Z4L  =  “ _  7 


3t  4ttcj 


2§ 


(44) 


In  general,  however,  all  of  the  terms  of  Equation  (43)  contribute  to  the 
generation  of  8,  although  the  contribution  from  higher-order  terms  may  be 
small  in  some  cases. 

Under  certain  circumstances,  the  electron  inertia  and  the  Hall 
terms  from  Equation  (39)  may  be  small.  In  this  case, 


1 

e(l+c)n. 


VP 
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The  last  term  in  the  above  equation  is  an  effective  electric  field  in  the 
direction  of  the  fluid  pressure  gradient — similar  to  the  thermoelectric 
effect.  The  curl  of  this  last  term  generates  a  magnetic  induction  in  the 
direction  of  Vng  x  vT.  In  other  words,  crossed  gradients  in  electron 
number  density  and  fluid  temperature  extract  magnetic  field  energy  from 
internal  energy  of  the  fluid. 

The  energy  equation  for  a  nonrelativistic  conducting,  viscous 
fluid  may  be  written  as: 


where  e  is  the  specific  internal  energy,  a1J  the  viscous  stress  dyad,  and  < 
the  thermal  conductivity  of  the  fluid.  Inserting  the  electric  field,  E, 
from  Equation  (45)  into  the  above  yields: 


de  1  .2 
p3t  *  a  J 


(<VT) 


VP 


(47) 


The  last  term  in  Equation  (47)  is  the  source  of  energy  for  thermally 
generated  magnetic  fields. 
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SECTION  5 

THE  THERMAL  GENERATION  OF  MAGNETIC  FIELDS 

5.1  INTRODUCTION. 

Recent  experiments  performed  at  NRL  with  the  Pharos  II  and  the 
Pharos  III  design  involved  the  illuminating  of  thin  planar  targets  by 
hundred-Joule  laser  pulses  with  intensities  between  IQ12  and  1014  Watt/cm2 
that  endure  for  a  few  nanoseconds.  The  computer  code  MACH2  (Peterkin,  et 
al.,  1986,  Ref.  7),  designed  to  perform  two-dimensional  simulations  of 
plasmas  in  complex  geometries,  was  used  to  simulate  some  of  the  NRL  experi¬ 
ments  that  took  place  in  environments  of  Nitrogen  gas  at  pressures  of  a  few 
Torr.  In  a  previous  report  (Stell ingwerf ,  Peterkin,  Longmire,  1986,  Ref. 
8),  we  presented  results  of  simulations  of  the  laser-target-ambient  system 
with  MACH2  that  illustrated  the  time  development  of  the  bubble  that 
emanates  from  the  laser-irradiated  target  and  the  laser-induced  cylindrical 
expansion  that  leads  to  the  formation  of  the  aneurysm.  In  this  section,  we 
investigate  the  formation  and  development  of  magnetic  fields’ that  are 
generated  from  thermal  energy  near  the  laser  focal  soot. 

5.2  LASER  AND  TARGET  PARAMETERS. 

The  purpose  of  our  investigation  is  two-fold.  First,  we  wish  to 
learn  the  magnitude  as  'well  as  the  spatial  and  temporal  development  of 
laser-induced  magnetic  fields.  Second,  we  wish  to  learn  what  effect  these 
fields  may  have  on  the  dynamics  of  the  laser-target-ambient  system.  To 
accomplish  this,  we  performed  two  simulations  from  the  same  set  of  initial 
conditions.  In  the  first,  we  do  not  nodel  thermally-generated  magnetic 
fields  and  in  the  second,  we  do.  We  then  compare  the  state  of  the  two 
simulations  at  identical  times  to  learn  what  effect,  if  any,  the  thermally- 
generated  magnetic  fields  have  on  the  dynamics  of  the  system. 
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Laser  and  target  parameters  used  for  the  MACH2  simulations  dis¬ 
cussed  herein  are  similar  to  those  of  a  particular  experiment  performed  at 
NRL  in  late  1985  (shot  85-544).  The  computational  mesh  used  for  the 
present  sequence  of  MACH2  calculations  is  illustrated  in  Figure  25.  The 
vertical  line  at  the  left  of  this  figure  is  an  axis  of  symmetry;  rotation 
about  this  line  is  assumed.  The  target,  located  at  the  bottom  of  the  mesh, 
is  100  urn  thick,  initially  of  density  2.7  gm/cc  and  at  room  temperature 
(0.025  eV).  The  ambient  nitrogen  pressure  is  5  Torr  and  extends  5  mm  from 
the  target  surface.  For  this  investigation,  we  are  interested  in  early- 
time  development  (the  first  15  ns).  To  obtain  sufficient  resolution  with  a 
modest  number  of  computational  cells,  we  restrict  our  attention  to  the 
region  close  to  the  target  surface.  The  laser  spot  is  a  circle  of  diameter 
100  urn.  The  laser  beam,  centered  on  the  axis  of  symmetry,  propagates  from 
the  top  of  the  mesh  to  the  critical  surface--that  is  the  target  surface  at 
early  times.  The  full  opening  angle  of  the  converging  laser  beam  is  5.92 
degrees.  The  computational  mesh  is  36  x  36  cells  with  four  cells  across 
the  beam  and  4  cells  across  the  thickness  of  the  target.  The  boundaries  of 
the  computational  region  allow  flow  of  material  through  them,  and  are 
electrically  non-conducting. 

The  target  in  shot  35-544  was  made  from  aluminium,  and  the  amoi- 
ent  gas  was  nitrogen  at  an  initial  pressure  of  5  Torr  at  roam  temperature. 
The  target  was  a  flat  aluminium  disk  with  a  radius  of  750  urn  and  a  thick¬ 
ness  of  5.6  ym.  The  thickness  should  have  little  effect  on  the  blowoff 
emanatinq  from  the  laser  side  of  the  target,  so  to  avoid  an  unnecessarily 
small  time-steo,  we  model  a  100  urn  thick  target  in  the  simulation  of  this 
shot. 


In  these  simulations,  a  pulse  of  1.054  urn  wavelenqth  laser  light, 
with  an  essentially  Gaussian  temporal  shape,  illuminates  the  target  for 
8.168  ns.  The  peak  power  occurs  at  4.084  ns,  and  the  full  width  at  half 
maximum  is  6.0  ns.  Peak  power  is  2.29  1010  Watt,  and  a  total  of  109  J  of 
energy  is  delivered  to  the  target.  The  peak  intensity  is  6  IQ13  Watt/cm2. 
The  laser  is  focused  onto  a  100  mn  radius  circle  on  the  target.  All  of  the 
energy  that  is  not  absorbed  by  the  ambient  gas  in  the  laser  channel  is 
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Figure  25.  MACH2  calculation  mesh  for  models  sgmf2  and  sgmf4 
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absorbed  at  the  critical  surface  (the  surface  in  three-dimensional  space  at 
which  the  laser  light  first  encounters  material  with  density  above,  the 
critical  density--3.435  10‘3  gm/cc  for  aluminium  and  1.054  um  light). 

Laser  light  is  assumed  to  be  absorbed  via  inverse  bremsstrahlung.  The  cold 

ambient  gas  is  preionized  in  our  models,  presumably  by  X-rays  from  the 

target  induced  by  the  laser  precursor.  Our  preionization  model  falls  as 
l/rz  from  the  target  surface  with  five  percent  of  the  ambient  atoms  sinqly 
ionized  1  cm  from  the  target.  The  preionization  does  not  affect  the 
temperature  of  the  ambient  gas. 

Thermal ly-generated  magnetic  field  arise  from  a  curl  of  the 
pressure  gradient  divided  by  the  sum  of  the  electron  and  the  ion  number 

densities.  Because  these  lie  in  the  computational  plane  of  MACH2,  the  curl 

is  orthogonal  to  that  plane  and  the  resulting  magnetic  field  is  in  the 
azimuthal  direction. 

5.3  SIMULATION  RESULTS. 

The  diagrams  in  Figure  26  illustrate  the  first  simulation  (with¬ 
out  a  model  of  thermally-generated  fields)  10  ns  after  the  laser  pulse 
beqan.  The  pictures  of  electron  number  density  and  temoerature  clearly 
show  the  cylindrical  expansion  from  the  laser  heated  channel.  As  the 
temperature  in  the  channel  rises,  the  density  there  falls  below  its  ambient 
value  in  order  to  maintain  approximate  oressure  equilibrium.  The  picture 
of  velocity  vectors  shows  a  jet  of  material  proDagating  up  the  low-density 
channel  ahead  of  the  expanding  bubble  outside  the  channel.  (It  is  worth 
noting  that  simulations  without  laser  heating  of  the  channel  show  no  such 
oerturbation  of  the  otherwise  spherical  bubble  that  emanates  from  the 
laser-heated  target.) 

The  diagrams  in  Figure  27  illustrate  the  second  simulation  (with 
a  model  of  thermally-generated  fields)  10  ns  after  the  laser  pulse  began. 

As  with  the  first  simulation,  the  pictures  of  electron  number  density  and 
temperature  clearly  show  the  cylindrical  expansion  from  the  laser  heated 


65 


generated  magnetic  fields. 


i/a «» 


Fiqure  27.  MACH2  model  sgmf4  at  10  ns  with  a  model  of  thermally-generated 
magnetic  fields. 


channel.  In  addition,  a  contour  plot  of  rB  (essentially  current  stream¬ 
lines)  is  shown  in  Figure  27.  Toroidal  regions  of  field,  with  a  peak  value 
of  37  kG,  appear  concentrated  directly  behind  the  bubble.  This  toroidal 
region  is  oriented  approximately  60  degrees  from  the  laser  axis.  Poloidal 
current  loops,  with  a  peak  current  of  15  kA,  circulate  about  the  toroidal 
field  lines.  At  this  time  100  kA  is  flowing  in  the  region  between  the  axis 
and  a  location  2  nm  from  the  axis.  At  5  ns,  immediately  after  the  time  of 
peak  laser  power,  the  magnetic  induction  reaches  a  maximum  of  nearly  1  MG 
in  a  limited  region  (one  computational  cell)  near  the  laser  focal  spot.  At 
this  time,  100  kG  fields  envelope  an  extended  region  (a  few  computational 
cells)  centered  about  the  focal  spot. 

Despite  the  size  of  the  current  and  magnetic  field  at  early 
times,  the  state  of  the  simulation  at  10  ns  is  essentially  unaffected  by 
the  thermal  generation  of  magnetic  fields.  This  leads  us  to  conclude  that, 
although  rather  strong  fields  are  generated,  they  are  a  small  perturbation 
and  do  not  change  the  dynamics  of  the  laser-target  interaction  by  a  notice¬ 
able  amount.  Although  multi-ki logauss  fields  are  generated,  the  fraction 
of  laser  energy  that  is  converted  to  magnetic  energy  peaks  at  less  than  one 
percent  as  can  be  seen  in  Figure  23  if  one  recalls  that  the  laser  energy  is 
109  J. 


Note  also  from  Figure  28  that  the  magnetic  energy  peaks  at  the 
end  of  the  laser  pulse  (8.2  ns).  After  this  time  the  magnetic  energy, 
which  arose  from  plasma  thermal  enerqy  in  the  first  place,  returns  to 
plasma  energy  by  the  Joule  heating  mechanism.  We  conclude  from  this  that 
the  magnetic  field  is  generated  not  by  dynamics  in  the  expanding  bubble, 
but  rather  by  the  laser  interacting  with  plasma  near  the  critical  surface. 
The  non-zero  size  of  the  laser  spot  and  the  extent  of  the  target  create 
two-dimensional  effects  that  are  important  to  the  early-time  evolution  of 
the  expanding  olasma  and  that  create  the  crossed  gradients  of  temperature 
and  density  that  give  rise  to  magnetic  fields. 
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The  expanding  bubble  is  hot  (tens  of  eV)  and  has  a  relatively 
high  electrical  conductivity.  The  magnetic  field  is  therefore  transported 
with  the  expanding  plasma  so  field  tends  to  reside  within  the  bubble.  This 
view  is  supported  by  Figure  29  where  we  show  the  state  of  the  problem  at 
1  ns  intervals  for  the  first  10  ns  of  the  second  simulation.  The  plasma  is 
not  diffusive  and  the  magnetic  field  is  essentially  trapped  within  the 
bubble. 


Our  model  for  thermal  conduction  is  isotropic.  However,  it  is 
well  known  that  thermal  conduction  in  a  direction  orthogonal  to  a  strong 
magnetic  field  (“eTg»l)  is  suppressed  because  of  the  small  gyroradius  of 
the  electron.  Figures  30  and  31  illustrate  the  same  simulation  at  10  ns, 
but  now  with  the  anisotropic  effects  included  (model  sgmflO).  The  overall 
plasma  conf iguration  is  very  similar  to  that  seen  in  Figure  27.  The  only 
significant  difference  between  the  two  runs  is  a  strong  suppression  of  the 
thermal  flux  from  the  heated  expanding  shell  in  sgmflO,  except  in  the 
heated  channel.  This  has  some  effect  on  the  temperature  distribution  in 
the  shell  and  in  the  channel,  but  has  little  effect  on  the  dynamics  of  the 
expanding  shell  and  the  aneurysm. 
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Figure  29  rB^  plots  at  1  ns  intervals  during  laser  pulse  from  MACH2 
simulation  sginf4. 


Fiqure  11.  Thermal  flux  and  temoerature  at  10  ns  for  the  model  sqmflO 
including  the  effects  of  anisotropic  maqnetic  fields. 


SECTION  6 


MAGNETIC  BUBBLE  FORMATION  IN  EXPANDING  LASER  PLASMA 

As  discussed  in  previous  sections,  the  2  1/2  0  MHO  code  MACH2  has 
the  capability  of  modeling  an  expanding  plasma  in  the  presence  of  a  complex 
magnetic  field  configuration.  Large  toroidal  fields  generated  by  thermal 
gradients  have  been  discussed  in  Sections  4  and  5.  It  is  of  interest  to 
inquire  whether  such  fields  will  interact  with  an  ambient  magnetic  field, 
and  what  the  evolution  of  the  ambient  field  would  be  in  such  a  case.  In 
this  section  we  present  the  results  of  initial  tests  of  a  laser/target 
simulation  with  a  10  kG  ambient  magnetic  field,  a  value  accessible  to  the 
NRL  experimental  setup.  This  run  is  designated  "laser29". 

The  geometry  of  the  simulation  is  shown  in  Figure  32.  The  param¬ 
eters  for  this  test  computation  are  listed  in  Table  3. 


Table  3.  Laser29  parameters. 


Laser  Energy 

I  in 

J 

Ambient  Pressure 

5 

Torr 

Pulse  FWHM 

6 

ns 

Spot  Size 

200 

mi  crons 

Ambient  Field 

10 

kG 

Grid  Height 

0.5 

cm 

Grid  Width 

0.25 

cm 

The  field  configuration  used  here  is  unusual:  the  ambient  field 
is  parallel  to  the  laser  beam.  This  configuration  has  been  chosen  to 
preserve  the  cylindrical  symmetry  assumed  in  the  simulation,  and  is  a 
possible  experimental  configuration.  Although  the  bulk  of  the  plasma 
motion  is  along  the  field  lines  in  this  case,  there  remains  a  substantial 
cross-field  motion  as  the  spherical  blast  wave  expands.  In  addition,  this 
case  allows  a  study  of  the  behavior  of  the  field  near  the  plasma/ 
(conducting)  target  interface. 
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Figure  32.  Geometry  used  in  the  ambient  8-fieid  simulation.  Parameters 
are  given  in  Table  3. 
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The  results  of  this  run  are  shown  in  Figures  33  through  37.  The 
overall  plasma  configuration  and  plot  diagnostics  are  very  similar  to  the 
models  discussed  in  previous  sections,  and  detailed  descriptions  will  not 
be  repeated  here.  Each  of  these  figures  shows  the  situation  at  5  ns,  just 
following  the  peak  of  the  laser  pulse  (at  4  ns),  and  also  at  10  ns,  just 
after  the  laser  intensity  drops  to  zero  (at  8  ns). 

We  note  from  Figures  33  through  37  that  the  development  of  the 
shell  and  aneurysm  is  unaffected  by  the  ambient  field  in  this  case.  A  com¬ 
parison  with  the  models  in  Section  5  indicates  that  the  evolution  of  the 
self-generated  magnetic  field  is  also  unaffected  by  the  ambient  field.  The 
configuration  of  the  field  itself,  shown  in  Figure  36,  shows  the  field 
lines  swept  aside  by  the  expanding  plasma,  while  maintaining  their  con¬ 
nection  at  the  target.  The  enhanced  field  at  the  axis  near  the  target  is 
an  error,  and  does  not  occur  in  problems  that  are  run  with  the  most  recent 
version  of  MACH2.  The  magnitude  of  the  ambient  field  is  not  greatly 
changed  by  the  motion  of  the  plasma,  and  remains  small  compared  with  the 
size  of  the  sel f-generated  fields  within  the  cavity. 

This  test  computation  is  intended  to  be  a  test  of  the  code  for 
this  type  of  problem,  and  to  give  some  hint  of  what  to  expect.  In  a  future 
report  we  will  discuss  the  physics  of  the  expanding  magnetic  cavity,  and 
present  a  full  series  of  models  of  various  cases. 
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Figure  33.  Velocity  vector  plots  for  the  run  “laser29' 
Maximum  velocities  at  the  two  times  are  37( 


Fiqure  34.  Temperature  contour  plots  for  the  run  Mlaser29"  at  5  and  10  ns. 

The  channel  temperatures  at  the  two  times  are  74  eV  and  47  eV. 


35.  Electron  density  contour  plots  for  the  run  ”laser29”  at  5  and 
10  ns.  Channel  densities  at  the  two  times  are  2.1elR  and 
3.4el7. 


Fiqure  36.  Maqnetic  field  vector  plots  for  the  run  "laser29,‘  at  5  and 
10  ns.  Maximum  values  at  the  two  times  are  36  kG  and  40  kG. 


Figure  37.  Toroidal  roaqnetic  field  contour  plots  for  the  run  ”laser29“  at 
5  and  10  ns.  Maximum  values  at  the  two  times  are  470  kG  and 
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